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ABSTRACT 


A  research  investigation,  both  experimental  and  analytical,  into  the 
generation,  suppression  and  prediction  of  fan  noise  has  been  completed. 

The  technical  effort  of  the  program  was  divided  into  three  major  areas: 
theoretical  investigations,  acoustic  normalizing  techniques,  and  noise 
minimization  and  suppression  design  techniques.  The  three  basic  areas 
of  investigation  were  supported  by  an  extensive  accumulation  of  scale 
model  and  full  scale  fan  noise  measurements,  which  have  been  presented 
in  this  report.  Significant  results  have  been  obtained  in  each  of  the 
areas. 

The  theoretical  investigations  have  provided  insight  into  the  basic  noise 
generating  mechanism  of  fan  noise  commonly  called  fan  or  compressor 
'‘whine" .  Analysis  of  the  mechanism  has  included  the  classical  theories 
of  sound  propagation  as  well  as  aerodynamic  cascade  theories  and  small 
perturbation  concepts. 

Acoustic  normalizing  techniques  have  resulted  in  a  correlation  curve  of 
the  overall  sound  power  for  13  separate  fern  and  compressor  vehicles  with 
each  one  being  from  a  different  design  family.  The  curve  covers  designs 
from  scale  model  to  full  scale  and  is  useful  in  determining  scaling 
effects  as  well  as  providing  a  basic  guide  to  design  engineers. 

Investigation  of  noise  minimization  through  design  techniques  incorporates 
the  results  of  the  analytical  investigations  and  the  normalized  sound 
power  relations.  Results  indicate  the  effect  of  rotor-stator  spacing, 
blade  number,  RPM,  inlet  length  and  centerbody  design. 

Suppression  techniques  evaluated  Were  wave  cancellation,  deflection,  and 
absorptive  suppressor  designs.  Test  results  of  an  absorptive  type 
suppressor  show  up  to  10.5  db  maximum  suppression  as  tested  with  a  scale 
model  fan  configuration. 

The  noise  data  accumulated  throughout  the  program  are  presented  in  table 
form  for  the  sound  pressure  measurements  made  on  each  configuration. 
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Recent  developments  in  advanced  propulsion  system  designs  have  utilized 
fan  type  configurations,  such  as  the  General  Electric/Ryan  XV5A  air¬ 
craft,  Figure  1  .  The  noise  associated  vith  a  system  of  this  type  is 
generally  less  on  a  sound  pressure  level  (SPL)  basis  than  that  associated 
vith  turbojet  systems  of  the  same  thrust  class  due  to  the  lover  Jet 
velocities.  Hovever,  the  fan  noise  or  "whine"  being  generated  by  the 
rotor  at  the  blade  passing  frequency  is  a  very  distinct  sound  that  re¬ 
quires  special  consideration  since  standard  Jet  noise  suppression 
techniques  are  not  applicable.  This  fan  "whine"  is  the  same  noise 
generating  mechanism  as  that  associated  with  compressor  noise.  Some 
analysis  of  this  type  of  discrete  tone  noise  mechanism  has  been  done 
previously,  but  only  a  limited  amount  of  published  data  is  available  for 
evaluation  of  various  analytical  techniques.  This  effort  and  report 
provide  extensive  supporting  data  for  analytical  investigations,  past  and 
present,  and  a  normalizing  technique  on  a  sound  power  level  basis  (PWL) 
for  the  performance  of  all  the  various  design  fans  and  compressors 
tested.*  In  addition,  suppression  concepts  have  been  investigated,  and  a 
10.5-db  reduction  absorptive  type  suppressor  that  will  provide  a  guide 
for  future  absorptive  type  suppressor  applications  has  been  designed  and 
demonstrated. 

This  report  is  divided  into  sections  corresponding  to  three  major  areas: 
analytical  investigations,  acoustic  normalizing  techniques,  and  noise 
minimization  and  suppression  design  techniques.  A  fourth  section 
describes  the  vehicles  tested  and  presents  the  basic  data  accumulated 
during  the  program. 


Figure  1.  XV- BA  Aircraft  in  Hover  Mode 
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SECTION  ONE  -  ANALYTICAL  INVESTIGATIONS 


Analysis  of  discrete  tone  noise  from  fans  is  made  difficult  by  the 
many  nonlinearties  involved  in  the  physical  aspects  of  the  problem. 
Previous  investigations  (references  15  and  18)  obtained  qualitative 
results  by  using  the  basic  wave  equation  with  appropriate  boundary 
and  Initial  conditions  and  making  necessary  assumptions  such  as 
rigid  walls,  small  pressure  fluctuations  and  homogeneous  medium. 

These  results  defined  the  sound  field  both  in  the  fan  duct  upstream 
and  downstream  of  the  rotor  face  and  in  the  radiated  field  at  a 
distance  from  the  fan  inlet  or  exhaust.  The  supporting  data  for 
these  investigations,  however,  was  sparse  and  did  not  adequately 
verify  the  results. 

The  problem  of  defining  the  discrete  tone  noise  generating  mechanism 
at  the  fan  rotor  is  generally  divided  into  two  cases:  rotor  or 
rotational  noise  associated  with  the  pressure  field  around  the  rotor 
blade,  and  rotor  -  stator  interference  noise  resulting  from  the  inter¬ 
action  of  the  viscous  wake  from  a  stator  row  upstream  of  the  rotor 
with  the  rotor  blades  or  the  viscous  wake  from  the  rotor  blade 
interacting  with  a  downstream  stator  row.  In  all  cases,  the  discrete 
tone  fundamental  frequency  or  first  harmonic  is  defined  by  the  product 
of  the  rotor  speed,  RPS,  and  the  number  of  rotor  blades,  Nr.  The 
analysis  of  this  mechanism  using  the  classical  wave  equation  results 
in  the  definition  of  a  "cutoff"  frequency,  which  is  the  frequency  above 
which  the  sound  waves  propagate  from  the  rotor  face  and  below  which  the 
sound  is  attenuated  at  a  very  fast  rate.  This  cutoff  frequency  is 
defined  for  both  the  case  of  rotational  noise  and  the  case  of 
rotor  -  stator  interference  noise.  Appendix  I  is  an  analysis  similar 
to  that  made  by  other  investigators  which  shows  the  development  of  the 
cutoff  frequency  concept. 

In  the  case  of  rotational  noise,  the  cutoff  frequency  is  not  attained 
until  rotor  tip  speeds  approach  Mach  1.0;  thus,  for  the  various  scale 
model  vehicles  tested  during  this  program,  the  necessary  cutoff 
frequency  was  not  reached.  The  rate  of  attenuation  predicted  for 
the  sound  field  below  the  cutoff  frequency  is  in  the  range  of  55  db 
per  blade  space  (reference  18,  page  60),  which  would  mean  that  the 
resultant  sound  pressure  levels  at  distances  of  ten  feet  or  more 
from  the  fan  face  would  be  insignificant.  This  same  rate  of  attenuation 
was  predicted  using  flat  plate  cascade  theory  as  shown  in  Appendix  II 
and  Appendix  III  and  Figures  2  and  5»  Figure  2  shows  the  attenuation 
in  velocity  for  the  limiting  case  of  no-flow  velocity,  while  Figure  5 
shows  the  effect  due  to  Mach  number,  M,  and  flow  angle,  p.  The  test 
vehicles.  Section  Four,  did  not  substantiate  this  theory,  since  all 
the  vehicles  tested  with  a  lone  rotor  produced  significant  sound 
pressure  levels  at  distances  up  to  50  feet. 
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Figure  2.1.  Velocity  Field  of  a  Flat  Plate  Cascade. 
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Figure  2.3.  Velocity  Field  of  a  Flat  Plate  Cascade 


Upstream 
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Figure  2.4.  Velocity  Field  of  a  Flat  Plate  Cascade. 
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Sound  Pressure  Level  Decay  per  Blade  Space  (db) 
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Mach  Number 

Figure  3.  Sound  Pressure  Level  Decay  vs  Tip  Mach  Number 
for  Various  Flow  Angles. 
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The  case  of  rotor  -  stator  noise  requires  a  variation  in  the  cutoff 
frequency  to  account  for  the  different  pressure  pattern  at  the  rotor 
face  due  to  the  alignment  of  the  rotor  blades  and  stator  blades.  This 
is  shown  in  detail  in  Appendix  I.  The  conclusion  drawn  by  this 
analysis  and  the  previous  analyses  was  that  rotor  -  stator  noise  would 
be  substantially  more  significant  than  the  noise  generated  by  a 
lone  rotor.  Test  results  of  the  rotor  -  stator  fan  and  single-stage 
compressor  test  vehicles.  Section  Four,  Table  1,  with  and  without  a 
stator  row  are  shown  in  Figure  4.  These  test  results  show  that  the  rotor 
only  condition  is  still  significant  relative  to  the  rotor  -  stator  case 
and  that  a  considerable  variation  is  found  in  the  db  reduction  obtained 
by  removing  the  stator  row.  The  rotor  -  stator  fan  attained  only  a 
4-db  reduction,  while  the  single-stage  compressor  attained  a  22-db 
reduction. 

The  above  discussion  for  the  rotational  end  rotor  -  stator  noise  genera¬ 
ting  mechanism  indicates  that  further  investigation  is  required  to 
better  define  the  amount  of  interference  required  to  generate  the  high 
noise  levels.  For  example,  discontinuities  on  the  fan  inner  surface  or 
upstream  and  downstream  support  struts  may  be  sufficient  to  interact 
with  the  rotor  to  produce  discrete  tones.  Tests  with  no-lift  devices, 
such  as  wires,  placed  near  the  rotor  face  would  be  an  example  of  the 
type  of  investigation  required. 

Appendix  IV  is  an  analysis  based  on  the  interaction  of  cylindrical 
waves  to  define  the  sound  field  resulting  from  rotor  -  stator  inter¬ 
ference.  The  particular  advantage  of  this  type  of  analysis  is  the 
elimination  of  the  cutoff  frequency  as  it  is  presently  defined  and  the 
relatively  simple  geometry  .of  the  analytical  model.  The  cylindrical 
waves  produced  by  each  rotor  -  stator  interaction  are  space  and  time 
sequenced  so  as  to  combine  to  produce  plane  wave  fronts  that  propagate 
without  attenuation.  The  disadvantage  of  this  approach  is  the  lack  of 
a  simple  acoustic  potential  as  is  associated  with  plane  or  spherical 
waves,  thus  requiring  the  application  of  the  unsteady  flow  method  of 
characteristics.  The  results  of  the  analysis  were  programmed  for 
computer  use  as  presented  in  Appendix  IV,  but  further  investigation 
is  required  to  obtain  realistic  numerical  results. 

Having  defined  the  noise  mechanism  at  the  rotor  face  and  the  prop¬ 
agation  characteristic  to  the  fan  inlet  or  exhaust,  the  radiation 
pattern  to  space  or  directivity  must  be  defined.  This  analysis  is 
presented  in  Appendix  V.  The  directivity  patterns  predicted  by  the 
analysis  are  a  function  of  the  number  of  rotor  and  stator  blades,  Nr 
and  Ns  respectively,  RPM,  and  physical  size  of  the  radiating  face.  The 
analysis  was  programmed  to  utilize  a  computer  in  making  the  necessary 
calculations.  Typical  results  from  the  program  are  shown  in  Figures  5»1 
through  Note  that  the  program  predicts  a  considerable  decrease  in 
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Rotor 

Axis 


No.  of  Rotor  Blades  (Nr)  ■  32 
No.  of  Stator  Blades  (Ns)  ■  36 
Rotor  Speed  m  1890  RPM  (Ns)  •  36 


Figure  5.1.  Predicted  Directivity  Pattern 
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Sound  Pressure  Level  Relative  to  Maximum  Level  (db) 


Rotor  No.  of  Rotor  Blades  (Mr)  «  32 

Axis  No.  of  Stator  Blades  (Ns)  m  32 

Rotor  Speed  ■  1000  RPM 


Figure  5.2.  Predicted  Directivity  Pattern. 
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Rotor 

Axis  No.  of  Rotor  Blades  (Nr)  *  32 


No.  of  Stator  Blades  (Ns)  «  33 
Rotor  Speed  ■>  1000  RPII 


Figure  5.3*  Predicted  Directivity  Pattern. 
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Sound  Pressure  Level  Relative  to  Maximum  Level  (db) 


Rotor  No.  of  Rotor  Blades  (Nr)  «  32 

Axis  No.  of  Stator  Blades  (Ns)  m  34 

Rotor  Speed  »  1000  RPM 


Figure  5.4  .  Predicted  Directivity  Pattern 
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sound  pressure  level  at  sane  positions.  This  is  the  sane  type  of 
directivity  predicted  in  reference  18,  page  35  .  Variations  in  sound 
pressure  level  in  the  order  of  15  db  are  not  readily  measured  in 
actual  cases,  as  evidenced  by  the  data  in  Section  Four  and  as  shown  in 
Figures  6  and  7.  Further  investigation  of  measured  directivity  patterns 
shoved  co-  ; 1  r  ble  variation  in  the  patterns  as  a  function  of  the 
distance  fan  face.  Typical  results  are  shown  in  Figures  8  and 

9.  The  radial  variations  do  not  correspond  to  the  spherical  divergence 
relationship,  as  3hovn  in  Figure  9*  lu  addition  to  the  variations 
found  in  each  test,  considerable  variation  is  found  in  measurements 
made  on  the  same  configuration  but  on  different  days.  Figure  10  showB 
the  variation  for  two  runs  made  on  different  days  with  the  development 
vehicle  fan,  Section  k,  Table  1.  These  inconsistencies  in  measuring 
directivity  patterns  make  correlation  of  analysis  and  test  results 
unrealistic.  In  addition,  they  raise  the  question  of  the  validity  of 
attempting  to  predict  directivity,  since  in  practical  applications 
the  many  variables  encountered  would  make  predictions  unrealistic. 


Figure  6.  Sound  Pressure  Level  Circumferential  Survey  of  CJ805-23  Fan, 
200-foot  Radius  (reference  Table  1) . 
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Figure  7  .  Rotor-Stator  Fan  Circumferential  Survey  (reference  Table  l). 


Figure  8.  Radial  Directivity  Variation  of  Development  Vehicle 
(reference  Table  1). 
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Sound  Pressure  Level  Circumferential  Survey  of  Development  Vehicle 
(reference  Table  1) . 


SECTION  TWO  -  ACOUSTIC  NORMALIZING  TECHNIQUES 


When  a  new  fan  design  is  in  the  development  stage  or  when  scale  model 
fan  test  results  are  being  extrapolated  to  predict  full  scale  per¬ 
formance,  applicable  normalizing  techniques  are  required.  One  of 
the  methods  used  previously  to  determine  acoustical  performance  of 
fans  or  compressors  was  a  curve  normalized  on  a  shaft  power  basis  and 
applicable  to  a  particular  design  family.  A  typical  prediction  curve 
for  various  families  is  shown  in  Figure  11.  This  curve  does  not  allow 
for  new  design  analysis  unless  a  similarity  is  maintained  between 
the  new  design  and  one  of  the  earlier  families  of  designs.  In  addition, 
little  or  no  data  was  available  to  substantiate  the  fact  that  tne 
families  of  design  lines  were  consistent  from  scale  model  to  full  scale. 
Utilizing  the  scale  model  and  full  scale  data  as  presented  in  Section 
Four,  shaft  power  was  used  to  normalize  the  fan  sound  power,  PWL,  for 
the  family  of  Vertical  Takeoff  and  Land,  VTOL,  vehicles :  inlet  guide 
vane  (IGV)  -  rotor  fan,  rotor  -  stator  fan,  lift  fan,  pitch  fan,  and 
cruise  fan,  all  of  which  are  described  in  Section  4,  Table  1.  This  in¬ 
dicated,  as  shown  in  Figure  12,  that  the  shaft  power  correlation  did 
hold  from  scale  model  to  full  scale  for  predicting  fan  or  compressor 
sound  power,  but  the  normalized  relationship  still  did  not  allow  eval¬ 
uation  of  new  designs. 

The  problem  of  normalizing  the  various  design  families  onto  one  curve 
was  solved  by  using  an  energy  flux,  E,  concept  and  a  design  parameter 

10  log  Aa  n  (Djj/Dj)2 


with  2 

A  =  Rotor  annulus  area,  ft 
a  7 

n  =  Rotor  speed,  RPM 

N  =  Number  of  rotor  blades 
r 

Dh/Dt  =  ratio  of  rotor  hub  diameter  to  rotor  tip  diameter. 

The  normalized  acoustic  performance  curve  is  shown  in  Figure  1.5 .  This 
empirical  curve  is  based  on  scale  model  and  full  scale  data  and  includes 
designs  ranging  from  the  62-inch  VTOL  lift  fan  to  scale  model  com¬ 
pressors.  The  solid  symbols  in  Figure  15  represent  data  accumulated  by 
the  program,  while  the  open  symbols  are  measurements  made  prior  to  the 
program.  The  abscissa  of  the  curve  represents  the  total  energy,  E,  of 
the  air  as  it  leaves  the  fan  or  compressor  rotor  stage  on  a  per-unit 
time  and  area  basis.  The  ordinate  of  the  curve  was  determined  by  a 
parametric  study  of  the  various  design  parameters  and  the  results  of  the 
noise  data  obtained  on  each  vehicle. 
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Mechanical  Power  of  Various  VTOL  Fans  vs  Sound 


Figure  13.  Normalized.  Overall  Power  of  Compressor  and  Fan  Noise 


Since  the  curve  is  based  on  the  fundamental  acoustic  parameter  sound 
power,  it  allows  the  designs  of  various  vehicles  to  be  compared 
directly.  This  type  of  analogy  is  useful  from  both  a  research  view¬ 
point  and  a  design  viewpoint.  For  research,  the  normalized  curve 
eliminates  many  of  the  irregularities  presently  found  in  fan  and  com¬ 
pressor  noise  measurements.  For  the  designer,  the  normalized  curve 
provides  a  basis  on  which  the  various  design  parameters  (rotor  annulus 
area,  Aa,  rotor  speed,  n,  rotor  blade  number,  Nj.,  hub  -  tip  ratio, 
Dg/Dp,  fan  air  flow,  W,  and  discharge  total  temperature,  Ty)  may  be 
evaluated  to  determine  the  optimum  combination  for  minimum  noise  gener¬ 
ation.  This  is  discussed  further  in  Section  Three. 

The  evaluation  of  advanced  designs  may  be  extended  from  the  sound 
power  level,  determined  by  the  normalized  power  curve,  to  a  sound 
pressure  level,  SPL,  by  using  additional  normalized  or  average  results 
from  the  test  data.  This  is  particularly  important  when  the  advanced 
design  must  conform  to  an  SPL  far  field  acoustic  requirement.  The 
overall  power  level  may  be  reduced  to  a  spectrum  power  level,  SPWL, 
for  each  harmonic  using  Figure  14.  This  normalized  power  spectrum  is 
the  result  of  averaging  the  measurements  taken  on  the  various  fan  con¬ 
figurations,  as  shown  in  Figure  15.  Note  from  Figure  15  that  a  con¬ 
siderable  spread  is  found  in  the  harmonic  power  spectrum  data.  Sound 
pressure  level  and  power  level  are  related  to  the  directivity  index, 

DI,  and  distance  from  the  source,  r,  by  the  following  equation: 

SPL  =  FWL  +  DI  -  20  log  r  -  10.5 

where  the  value  of  r  is  the  radius  at  which  the  SPL  is  to  be  determined 
and  the  value  of  the  directivity  index,  DI,  is  assumed  to  be  5*  The 
value  5  corresponds  to  the  average  DI  at  the  angle  of  maximum  noise  for 
the  various  vehicles  tested.  The  angle  of  maximum  noise  or  directivity 
was  not  normalized;  thus,  the  SPL  calculated  can  only  be  assumed  to  be 
in  the  vicinity  of  30  to  60  degrees  from  the  inlet  or  exhaust  of  the 
vehicle.  This  range  for  the  angle  of  maximum  noise  is  a  reasonable 
value  based  on  the  data  given  in  Section  Four.  If  a  design  is  to  be 
considered  that  is  similar  to  one  of  the  test  vehicle  configurations, 
a  more  realistic  value  for  the  angle  of  maximum  noise  may  be  determined 
from  Section  Four. 


As  an  illustration  of  the  procedure  discussed,  assume  the  following  fan 
design  parameters: 


Outer  diameter 

Inner  diameter 

Weight  flow 

Stage  temperature 
rise 


H 


W 


40  Inches 
30  inches 
150  lb/ sec 


AT  =  15°R 
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Spectrum  Sound  Power  Level  Relative  to  Overall  Sound  Power  Level 

SPWL  -  PWL  (db) 


1  2  3  4  5 

Harmonic  Number 


Figure  15.  Accumulated  Spectrum  Sound  Power  Level  Data 


RFM  n  =  o,000  RPM 

Rotor  blade  no.  I>  =  5^ 

Inlet  temperature  T  =  520 °R 

Calculate  the  following: 


dh/dt  =  30/40  =  .75 


(D^)2  =  ( .75)2  =  .562 


A  = 


it 

T 


x  (Dj‘ 


,3 

dt 


■T*©  L1  --562] 


“t 


A  =  -f-  x  11.2  x  .458  =  3.86  ft 
a  4 

Tt  =  T  +  AT  =  520 °R  +  15°  =  535 °R 
Total  Enthalpy  =  12b  at  535 °R  (from  gas  tables) 


H  x  W  „ 

-  ^r—  -  -  4.98  x  io3  7 


Btu 
sec-f 


From  Figure  13,  at  E  =  4.96  x  10^ 


Btu 


sec-ft2 


PWL  -  10  log 


W3sY 

N.  \DT/ 


2 

=  133  db 


10  log  (D^D  f  =  10  log  3  —  ^r-000  x  .562  =  10  log  322 

T  =  10  X  2.508  =  25 

PWL  =  133  +  25  =  158  db  overall 
Spectrum  PWL,  Figure  14: 

First  harmonic  =  PWL  -  3-5  =  lp8  -  3*5  =  15^.5  db 
Second  harmonic  =  PWL  -  5.5  '  15&  -  5*5  =  152.5  db 

Sound  pressure  level  at  angle  of  maximum  noise  (assume  a  100-foot  radius) 

SPL  =  PWL  +  DI  -  20  log  r  -  10. 5 
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SPL  =  158  ♦  5  -  20  log  100  -  10.5 
=  165  -  40  -  10.5  =  lfc5  -  50.5 
SPL  =  112.5  db  at  100  feet,  angle  of  maximum  noise 

The  above  procedure  allows  &  direct  analysis  of  the  acoustic  performance 
of  a  development  vehicle  base!  on  the  fan  design  parameters.  The 
question  of  installation  effects,  however,  requires  further  analysis. 
Several  installation  parameters  were  investigated  empirically:  Inlet 
length,  exhaust  duct  length,  centerbody  shape,  and  exhaust  louver 
vectoring. 

The  tests  were  made  by  varying  the  various  parameters  on  the  26- inch- 
diameter  VTOL  scale  model  fans  (reference  Section  Four,  Table  l). 

Figures  16  through  19  show  the  results  of  these  tests.  Hie  two  center- 
body  configurations.  Figure  16,  showed  no  effect  due  to  the  throttling 
of  the  longer  circular  centerbody  on  the  sound  power  distribution.  If 
the  circular  centerbody  was  sufficiently  large  in  diameter  to  block  the 
flow,  obviously  a  definite  effect  would  be  obtained.  Three  inlet 
lengths  were  tested:  4  inches,  9  inches,  and  23.5  inches.  Some  data 
scatter  is  evident  in  Figure  17,  but  no  definite  trends  are  indicated  to 
show  an  effect  due  to  the  varying  inlet  length.  Figure  18  shows  a  com¬ 
parison  of  two  configurations  with  and  without  a  l6-inch  exit  duct.  The 
two  configurations  also  vary  in  inlet  length  and  centerbody,  but  the 
previous  discussion  allows  a  direct  comparison.  The  data  as  shown  on 
Figure  18  has  some  scatter,  but  no  effect  due  to  the  long  exit  ducts 
may  be  concluded.  Results  of  varying  the  exhaust  louver  angle  with 
respect  to  the  flow  direction  also  showed  data  scatter.  Figure  19,  hut 
again  no  definite  trends  were  indicated.  Thus,  of  all  the  parameters 
varied, it  was  generally  concluded  that  no  definite  effect  could  be 
associated  with  installation  characteristics  excluding  radical  changes 
or  those  that  varied  the  fan  performance. 
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Figure  16.  Spectrum  Sound  Power  Level  Variation  with  Exhaust  Plug  Design 


Rotor-Stator  Fan  (reference  Table  1) 


VTOL  IGV  -  Rotor  Fan  (reference  Table  1) 
0°  IGV  Setting,  9-Inch  Inlet 
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Figure  19.  Spectrum  Sound  Power  Level  Variation  with  Exhaust  Louver  Vectoring  Angle 


SECTION  THREE  -  NOISE  MINIMIZATION  AND  SUPPRESSION 


MINIMUM  NOISE  GENERATION  PARAMETERS 
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Variovis  design  parameters  are  of  particular  interest  when  defining  a 
new  fan  configuration  for  minimum  noise  generation.  A  basic  design 
guide  to  be  followed,  at  this  point  in  the  state  of  the  art,  is  given 
by  Figure  13,  which  was  discussed  in  Section  Two.  This  normalized 
power  level  curve  shows  that  a  large  fan  doing  the  same  work  as  a 
smaller  fan  should  be  inherently  quieter  on  an  energy  flux  basis,  but  a 
small  fan  doing  proportionally  the  same  work  as  a  large  fan  may  generate 
nearly  the  same  amount  of  noise.  For  fans  of  a  given  power  requirement, 
the  parameter 

10  log  Aa  n  (Dg/Aj,)2 


may  be  evaluated  to  determine  qualitatively  the  proper  design  approach. 
Factors^in  the  numerator  of  the  normalizing  function,  A  ,  n,  and 
(Du/El,)  ,  should  be  minimized.  The  area,  A&,  also  appears  in  the  energy 
flux  term,  the  abscissa  of  the  curve.  Thusf  if  the  area  is  varied,  the 
amount  of  work  done  by  the  fan  must  be  varied  to  maintain  constant 
energy  flux;  but  this  requires  that  the  fan  performance  requirement  be 
redefined.  The  factor  n,  or  RPM,  shows  that  the  work,  if  done  at  a 
slower  speed,  will  generate  less2noise,  which  is  consistent  with  pro- 
pellor  theory.  The  term  (D„/lL,)  ,  hub-tip  ratio  squared,  is  a  measure 
of  the  radial  distribution  of  the  blade  loading.  The  higher  D^/D 
berms  are  representative  of  a  short  blade;  thus,  for  the  same  number  of 
blades  and  the  same  work,  they  have  a  higher  blade  loading.  With  the 
lower  D.  /D  terms  is  associated  a  smaller  blade  loading  on  the  same 
basis..  The  denominator,  N  ,  blade  number,  has  the  same  correlation  as 
Djj/D  ,  in  that  the  higher  £lade  number  for  the  same  fan  power  indicates 
a  lower  blade  loading  and  obviously  the  smaller  blade  number  indicates 
a  higher  blade  loading.  This,  too,  is  consistent  with  propellor 
theory.  Utilization  of  the  normalized  power  curve  in  the  above  manner 
could  provide  at  least  a  guide  in  defining  the  most  optimum  configura¬ 
tion  acoustically.  In  particular,  designs  that  are  near  the  knee  of 
the  curve,  energy  flux  values  up  to 

3  x  103  Btu 

2 

sec  -  ft 

could  benefit  appreciably  from  proper  design. 
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One  of  the  cannon  characteristics  of  all  the  various  fan  and  compressor 
designs  shown  on  Figure  13  is  the  relatively  close  rotor  -  stator 
spacing.  This  spacing  has  a  definite  effect  on  the  noise  generated  by 
the  fan,  as  discussed  in  Section  One.  Several  tests  were  made  on  the 
scale  model  fan  configurations  to  define  the  amount  of  noise  reduction 
obtained  by  increasing  the  rotor  -  stator  gap.  Obviously  the  case  of  a 
lone  rotor  is  the  maximum  amount  of  noise  reduction  that  can  be  obtained 
through  this  particular  design  parameter  variation.  Figure  20  shows 
the  results  of  these  tests  in  terms  of  power  level  variation  with  social 
rotor  -  stator  spacing.  The  extreme  variation  in  maximum  reduction 
obtained  without  a  stator  row  indicates  that  more  data  must  be  obtained 
to  determine  how  much  interference  is  required  to  produce  the  high 
sound  pressure  levels.  For  instance,  even  with  the  stators  removed, 
there  may  be  a  disturbance  between  the  rotor  and  some  discontinuity  in 
the  surface  of  the  stator  row  or  perhaps  a  disturbance  between  a  down¬ 
stream  or  upstream  strut  that  is  significant. 


Figure  20.  Variation  in  Sound  Power  Level  with  Rotor-Stator  Axial 
Spacing. 


38 


A  design  parameter  of  considerable  Interest  is  the  blade  passing 
frequency  of  the  fan,  the  product  of  the  number  of  rotor  blades  and  the 
RFM.  If  this  frequency  is  put  beyond  the  audible  range  through  high 
blade  numbers  and/or  high  RFM,  obviously  the  result  would  be  a  quiet 
fan.  Investigation  of  this  approach  was  beyond  the  scope  of  the  research 
investigation,  but  should  be  kept  in  mind  as  a  possible  means  of  elim¬ 
inating  noise. 


NOISE  SUPPRESSION  SCHEMES 


Three  suppression  schemes  were  evaluated  on  an  engineering  assessment 
basis:  deflection  methods,  wave  cancellation,  and  absorptive  suppression. 
Of  the  three  methods  assessed,  one,  absorptive  suppression,  was  selected 
for  further  evaluation  and  was  demonstrated  on  a  scale  model  fan  config¬ 
uration. 


DEFLECTION  METHODS 


The  concept  of  altering  the  directivity  pattern  of  a  particular  fan 
configuration  to  minimize  noise  in  specified  areas  was  assessed  qualita¬ 
tively  as  being  unrealistic  under  actual  operating  conditions. 

The  tests  described  in  Section  Two,  in  which  various  installation  charac¬ 
teristics  were  investigated,  showed  no  consistent  effect  on  the  directiv¬ 
ity  patterns,  while  results  of  Section  One  point  out  the  irregularity  of 
directivity  patterns  for  the  same  design  and  the  difficulty  in  predic¬ 
tion  of  directivity.  The  control  of  these  directivity  patterns  would 
require  a  major  effort  in  acoustical  design  for  each  operating  condition 
and  fan  configuration.  As  additional  information  is  obtained  with 
regard  to  the  basic  noise  generating  mechanism  of  fans  and  as  an  accu¬ 
rate  representation  of  the  sound  field  being  produced  by  the  rotor  is 
defined,  control  or  prediction  of  the  directivity  patterns  will  become 
more  realistic. 


WAVE  CANCELLATION 


The  pure-tone  characteristics  of  the  fan  noise  allows  an  analysis  to  be 
made,  assuming  a  sinusoidal  wave  shape,  that  shows  considerable  reduc¬ 
tion  in  level  obtained  by  varying  the  phase  of  two  signals  being  radi¬ 
ated  from  adjacent  sources.  Figure  21  shows  the  relative  amplitude  of 
the  addition  of  two  signals  as  a  function  of  the  phase  shift  between 
them.  Obviously,  a  180-degree  phase  shift  produces  total  cancellation; 
however,  with  a  shift  of  approximately  115  degrees,  a  reduction  in 
amplitude  of  50  percent  is  obtained  that  would  provide  a  notable  reduc¬ 
tion  in  noise  level.  The  'phase  shift  obtained  between  two  signals 
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Figure  21.  Relative  Amplitude  vs.  Angle  for  Various  Phase  Shifts. 

Variation  in  Wave  Amplitude  vs  Phase  Shift  for  the 
Addition  of  Two  Sine  Waves. 


aay  be  analyzed  as  a  function  of  the  distance  ,  L,  traveled  by  the  two 
waves  to  the  point  at  which  they  interact.  Figure  22  is  a  resultant 
plot  of  this  type  of  analysis  showing  the  difference  in  distance,  A  L, 
versus  frequency  for  various  degrees  of  phase  shift. 


fr^j< >c jf  a  lo’5  fCyclv*  pm r  m>coo4) 


Figure  22.  Correlation  of  Length  of  Wave  Travel,  Frequency,  and  Phase 
Shift  for  Sine  Wave  Interaction. 

Tests  were  conducted  on  a  system  of  this  type,  as  shown  in  Figures  23 
and  24.  Microphones  1  and  2  were  used  to  deter nine  phase  shifty  while 
microphone  3  was  used  to  determine  the  db  reduction.  The  results  are 
shown  in  Figure  25  in  terms  of  the  A  db  obtained  and  phase  shift  as 
compared  to  the  theoretical  reduction.  The  data  were  consistent  and 
showed  that  reduction  could  be  obtained  over  the  local  area  of  inter¬ 
ference  between  the  microphones;  however,  measurements  made  at  various 
distances  from  10  feet  to  1  foot  away  from  the  tubes  did  not  show  re¬ 
duction.  In  these  laboratory  tests,  the  phase  shift  was  obtained  by 
varying  the  distance  of  the  source  from  the  plane  of  emission,  but  in 
an  actual  application  this  is  not  practical;  thus,  the  phase  shift 
must  be  obtained  aerodynamically.  This  may  be  done  by  designing  al¬ 
ternate  diffusing  and  converging  passages  for  the  flow,  thus  varying 
the  velocity  at  which  the  signal  passes  through  each  section.  At  the 
same  frequency  and  a  different  velocity,  the  wave  length  varies;  thus, 
a  phase  shift  is  obtained.  The  system  used  for  demonstration  is 
shown  in  Figures  26  and  27.  This  flow  passage  was  designed  for  a  180- 
degree  phase  shift  at  an  inlet  Mach  number  of  .30  and  a  source  of 
.1670  cycle  per  second,  c/s. 
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These  figures  were  based  on  the  following  relationships  for  the 
two  pos sages : 


and 


at  l80  degrees  phase  shift 
L 


where 

f  =  frequency,  c/s 

At  =  time  difference,  sec 

t  =  time,  sec 

L  =  length  of  travel,  ft 


=  average  velocity  through  diffusing  passage,  ft/sec 
Vg  =  average  velocity  through  converging  passage,  ft/sec 


No  substantial  data  was  obtained  with  this  system  due  to  instability  of 
the  air  flow  through  the  passages,  the  difficulties  in  maintaining  a 
signal  that  could  be  effectively  changed  in  wave  length,  and  reflections 
of  the  signal  from  the  walls  of  the  flow  passage  sufficient  to  affect 
the  wave  length  variations. 


Ihe  conclusion  drawn  from  the  tests  was  that  the  wave  cancellation 
technique  would  require  excessive  aerodynamic  design  and  alteration  of 
the  fan  design,  such  as  "egg- crating"  the  exhaust  to  provide  significant 
noise  reduction. 


ABSORPTIVE  SUPPRESSION 


The  use  of  acoustical  material  to  suppress  fan  whine  had  been  demon¬ 
strated  previously  in  the  development  of  an  absorptive  suppressor  for 
the  CJ805-23  aft  fan  engine,  reference  Table  1.  This  suppressor  con¬ 
sisted  of  a  fiberglass-type  material  covered  with  a  fine  mesh  screen 
and  perforated  metal  and  located  on  the  inner  diameter  of  the  fan  exhaust 
passage.  The  suppression  obtained  was  quite  significant,  particularly 
on  a  A  db  per  linear  length  of  suppressor  basis,  being  approximately  1 
db  suppression  per  inch.  This  type  of  design  is  applicable  to  inlet  and 
outlet  guide  vane  modification  as  well  as  flow  passage  applications, 
assuming  that  this  same  suppression  rate  .nay  be  maintained  or  improved. 
This  method  of  suppression  was  chosen  for  further  investigation  and  for 
demonstration  on  a  scale  model  fan  configuration. 
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Figure  27.1  Flow  Ducts  for  Wave  Cancellation  Evaluation 
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ABSORPTIVE  SUPPRESSOR  DESIGN 


Absorptive  suppression  lends  itself  to  several  areas  of  physical 
application:  outlet  guide  vanes,  OGV,  inlet  guide  vanes,  IGV,  flow 
passages,  and,  as  in  the  case  o*'  an  "in-wing"  lift  fan  device,  exhaust 
vectoring  louvers.  This  latter  area  was  chosen  for  design  of  an 
absorptive  suppressor,  as  it  could  have  a  direct  application  to  develop¬ 
ment  vehicles  and  was  readily  adaptable  to  the  VTCL  scale  model  rotor- 
stator  fan  with  exhaust  louvers.  Figure  28. 


Figure  28. 


Exhaust  Vectoring  Louvers. 


DESIGN  METHOD 


Absorptive  suppressors  can  be  designed  readily  by  using  a  procedure  as 
outlined  in  reference  9>  page  243.  The  procedure  utilizes  a  design 
curve  shown  schematically  in  Figure  29.  The  design  parameters  are 
passage  width,  1*,  thickness  of  material,  t,  percent  open  area  (which 
are.  interrelated},  and  frequency.  By  varying  these  parameters,  a  specific 
flow  resistance,  R  ,  can  be  determined  for  the  thickness  of  material; 
and  from  this  dataxthe  proper  material  can  be  selected  by  using  a  curve 
similar  to  Figure  30.  This  general  outline  of  the  procedure  gives  some 
insight  into  the  problem  of  applying  this  type  of  design  to  a  passage 
where  there  is  air  flow.  If  the  aerodynamic  performance  is  to  be 
maintained,  the  proper  flow  conditions  must  be  kept;  but  they  are  in 
general  a  function  of  the  passage  width,  lengtl^  and  percent  open  area. 
ThuSy  a  compromise  between  aerodynamic  and  acoustical  design  must  be 
made. 


Figure  29.  Design  Chart  for  Acoustically  Treated 
Ducts  Lined  on  Two  Sides. 
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Lines  A  through 
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For  fans  and  compressors,  a  basic  design  parameter  for  stator  blades 
is  (reference  Figure  3l)* 


cr _  i 

AR  cos  ~  1 


Zs 


c/s 


h/c  cos  P 


where 


0O  =  flow  angle  out  of  blades 


2 

c 


m 


chord  length 

flow  passage  area  between  blades  at  stage  exit 


solidity  =  Wh  ■ 

J  blade  spacing 


£ 

S 


, .  blade  height  h 

aspect  ratio  =  — — -r—z — “rr:  =  —  • 

*  chord  length  c 


/$  -  exit  flow  angle 


Figure  31.  Schematic  of  Compressor  Blade  Design  Parameters. 
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Obviously,  as  the  blade  thickness,  t,  is  varied  for  addition  of  absorptive 
material,  the  chord,  c,  must  also  be  varied  to  maintain  the  thickness 
ratio,  t/c,  and  to  provide  sufficient  length  for  ample  suppression. 

With  a  varying  chord  length,  the  above  parameter  must  also  be  adjusted. 
Generally,  c  is  increased;  thus,  since  blade  height,  h,  is  predetermined, 
the  spacing, s,  must  be  increased  to  maintain  c/hs  a  constant.  When  s 
is  increased,  the  proper  percent  open  area, 

-  x  100, 
s 

must  be  maintained  to  satisfy  the  acoustical  design  parameters.  This 
type  of  design  procedure  becomes  difficult  when  a  low-frequency  signed 
is  being  attenuated  in  a  snail  passage,  since  the  material  thickness 
is  in  general  large  for  low-frequency  suppression  and  the  small  open 
area  is  affected  to  a  large  degree  by  the  increased  blade  thickness. 

For  the  application  of  absorptive  material  to  the  exhaust  louvers  of 
the  scale  model  fan  configuration,  the  fundamental  frequency  was  too 
low  to  provide  a  realistic  louver  configuration  due  to  the  large  material 
thickness  required.  To  demonstrate  a  suppressor  properly,  the  desig- 
was  made  to  suppress  the  higher  harmonics  of  the  fan  which  would 
correspond  to  the  fundamental  frequency  of  a  full  scale  vehicle.  The  de¬ 
sign  as  made  on  this  basis  would  also  be  of  the  same  general  physical 
configuration  as  would  be  required  by  a  full  scale  vehicle.  The  ab¬ 
sorptive  louver  design  is  shown  in  Figure  32,  and  as  mounted  on  the 
scale  model  fan  in  Figure  35.  The  design  calculations  are  given  in 
Appendix  VI. 

TEST  RESULTS 


The  absorptive  louvers  were  tested  Initially  in  a  reverberation  room 
environment  without  air  flow.  The  louvers  were  attached  to  a  large 
plenum.  Figure  33 >  by  using  a  loudspeaker  and  random  noise  generator 
as  a  sound  source,  the  signal  was  measured  for  the  standard  louver 
case.  The  tests  were  then  repeated  for  tne  treated  louver  case.  The 
results  are  shown  in  Figure  34  as  a  function  of  the  SPL  reduction 
obtained.  Also  shown  on  Figure  34  is  the  design  suppression.  The  re¬ 
sults  were  within  two  db  of  the  predicted  suppression  at  the  maximum 
suppression  point. 

The  louvers,  treated  and  standard,  were  next  tested  in  a  free  field 
environment,  Figure  35*  using  the  scale  model  rotor  -  stator  fan  to 
determine  the  effect  of  air  flow.  The  effect  of  air  flow  on  the 
suppression  characteristics  of  absorptive  suppressors  is  generally  to 
reduce  the  peak  suppression  by  25  percent  and  to  shift  the  point  of 
maximum  suppression  to  a  higher  frequency  (reference  17,  page  4).  The 
test  data.  Figure  34,  in  terms  of  pressure  level  reduction  indicated 
approximately  a  35  percent  reduction  in  suppression  and  a  shift  in  the 
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maximum  suppression  point  to  a  lower  frequency.  The  shift  to  a  lower 
frequency  was  not  expected  and  is  not  explainable  in  terms  of  the 
characteristics  of  absorptive  material.  The  results  show  a  very 
effective  reduction  in  noise  level  for  the  design,  and  the  maximum 
suppression  occurs  in  a  region  that  would  correspond  to  the  fundamental 
frequency  of  a  full  scale  vehicle.  The  variation  in  the  radiated  sound 
pressure  due  to  the  absorptive  suppressor  is  shown  in  Figures  30. 1 
through  36.*  for  the  first  5  harmonics.  As  would  be  expected,  th«. 
largest  suppression  occurs  on  the  exhaust  side  cf  the  fan. 
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Leading  Edge 


No.  Lorvers  I  X  Dimension, 
Required  for  |  in.  (Reference 
Figure  35) 


I 


Description 


Perforated  Metal,  20  Ga. 
1/ii"  holes  on  3/8"  centers. 


Flat  She°t,  l/32"  thick. 


,00Ji5  Wire  Screen,  80  mesh. 


Absorptive  Material , 
41b/ft3; 

R  =  1.5  x  10  n;ks  rayles 
meter 


Figure  32.  Absorptive  Louver  Design  for  Rotor-Stator 
Fan  (reference  Table  1)  . 


Figure  33.  Exhaust  Vectoring  Louvers  on  Reverberation 
Room  Plenum. 
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Figure  35.  Rotor>Stator  Fan  with  Treated  Louvers 
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Figure  36.1  Directivity  Pattern  of  Treated  and  Untreated  Louvers 
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Figure  36.2.  Directivity  Pattern  of  Treated  and  Untreated  Louvers. 
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Figure  36.3  Directivity  Pattern  of  Treated  and  Untreated  Louvers 


SECTION  FOUR  -  TEST  DATA 


TEST  VEHICLES 


The  following  data  describes  the  various  test  vehicles  used  during  the 
program  to  provide  information  as  described  in  previous  sections.  The 
data  given  is  for  the  basic  test  vehicle  and  does  not  account  for  various 
design  changes  as  they  are  presented  in  the  appropriate  sections.  The 
data  is  summarized  in  Table  1. 

VTOL  SCALE  MODEL  FANS 


Two  scale  model  fans  of  the  same  design  family  as  the  full  scale  VTOL 
lift  and  pitch  fans  used  in  the  XV5A  aircraft,  Figure  1,  were  tested. 
They  are  designated  the  IGV-rotor  fan  and  the  rotor  -  stator  fan,  the 
first  configuration  being  one  that  induces  swirl  ahead  of  the  rotor 
and  the  latter  being  one  that  induces  swirl  in  the  rotor. 

Inlet  guide  vane  -  rotor  fan  (Figure  37) 

D^,  =  26  inches 

DH 

r  *  *u5 

T 

W  =27-8  lb/sec 

RPM  =  1780  (constant  speed  drive) 

N  =  47 
r 

N  =  48 
s 

Pr  =  1.011  (fan  pressure  ratio) 
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TABLE  1 


UMMARY  OF  TEST  VEHICLE  CONFIGURATIONS 


Vehicle!  Tet ■ ed  During  the  Research  Investigation 


Name 

Inches 

BA 

w 

lb/aec 

RPM 

N 

r 

N 

s 

P 

r 

IGV -Rotor  Fan 

26 

.45 

27.8 

1780 

(Constant) 

47 

48 

1.011 

Rot or -Stator  Fan 

26 

.40 

2t  1 

1780 

(Constant ) 

36 

90 

1.005 

Lift  Fan 

62.5 

.40 

'2? 

2640 

(Maximum) 

36 

90 

1.115 

Pitch  Fan 

36 

■45 

144 

4074 

(>:-xloum) 

36 

90 

1.06 

Cruise  Fan 

62.5 

.40 

529 

264C 

(Maxima.') 

36 

90 

1.115 

Single  Stage  Scale 
Model  Compressor 

Development  Vehicle 

26 

.846  11 

Proprietary  Information 

1780 

(Constant) 

56 

48 

1.008 

Vehicles  Tested  Prior  to  the  Research  investigation 


Name 

Report  from  which  Data  was  Obtained 

Laboratory  Compressor 

Hunter.  W.J.  and  Smith,  L.H.,  Jr.,  Inlet  Guide  Vane-Rotor  Interference  Noise  Mens urement3, 
TIS  DF63FPD55,  General  Electric  Company,  Cincinnati,  Ohio,  January,  1963. 

CJ805-3  Compressor 

Semrau,  W.R.,  Effect  of  Variation  of  Stator  Vane  Angle  on  CJ805  Compressor  Noise,  DIM  425, 
General  Electric  Company,  Cincinnati,  Ohio,  May,  1961. 

J93  Compreaf.or 

Data  not  laiued  as  report. 

Research  Cosipressor 
No.  1 

Smith,  E.B.,  Noise  Measurements  of  MJ254  Compressor  FSCT  at  Lynn,  DIM  448,  General  Ele  trie 
Company,  Cincinnati,  Ohio,  November,  1961. 

Research  Compressor 
No.  2 

Lee,  R.  and  Semrau,  W.R.,  Low  Speed  and  Transonic  Speed  Research  Compressor  Acoustics 
Measurements,  DIM  457,  General  Electric  Company,  Cincinnati,  Ohio,  December,  1901. 

R.  Co.  12  Compressor 
and  RA26  Compressor 

Coles,  G.M.,  Title  Unknown,  Rolls  Royce  Company,  January,  I960. 

CJ805-23  Fan 

Dalke,  C.A.  and  Smith,  E.9. ,  Third  Noise  Measurement  of  the  GE  CJ805-23  Engine  and  Effects 
of  OGV's  and  Secondary  Nozzle  Sound  Treatment  on  Aft  Fan  Whine,  General  Electric  Company, 
Cincinnati,  Ohio,  November,  1959; 

Hallbach,  J.R.  and  Keener,  J.R.,  Basic  Design  Certification  Verification  Report  for 
7E-SSF-30  Sound  Suppressor,  TIS  R61FPD108,  General  Electric  Company,  Cincinnati,  Ohio, 
March,  1961. 

CF700  Fan 

McCann.  E.O.  and  Smith,  E.B.,  Noise  Measurement  of  the  CF700-2B  Engine,  DIM  534,  General 
Electric  Company,  Cincinnati,  Ohio,  June,  1963. 

Window  Fan 

Wells,  R.J.,  Industrial  Acoustics  Course  II  Lecture  No  es,  Fan  Nctse  Problem,  R53GL241-6, 
General  Electric  Company,  Schenectady,  N.Y.,  May,  1954. 
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Exhaust  Louver 
Vectoring  Angle 


Figure  37. 


Test  Conf igurations  of  Inlet  Guide  Vane  -  Rotor 
Fan  (reference  Table  1) . 


Rotor  -  stator  Fan  (Figures  3 8>  39  and  40) 


D_  =  26  inches 

T 

DH 

T 

W  =  26.1  lb/sec 

RFM  =  1780  (constant  speed  drive) 

N  =36 
r 

N  =90 
s 

P  =  1.005 

V 


Figure  38.  Rotor-Stator  Fan 
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23.5 
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stator  How 


3/8  Normal  Rotor-Stator 
Spacing 


Figure  40.  Test  Configurations  of  the  Rotor-Stator  Ean  (reference  Table  1) . 
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VTOL  LIFT  FAN 


The  VTOL  lift  fan  is  utilized  in  the  XV5A  aircraft  and  was  tested  as 
mounted  in  the  aircraft  and  in  a  static  test  facility  (Figures  1  and 

*L). 


fyj,  =  62.5  inches 

dh 

—  =  .40 

dt 

W  =  529  lb/sec  @  2640  RFM 

RPM  =  2640  maximum 

N  =36 
r 

N  =90 
s 

P  =  1.115 

r 

VTOL  PITCH  FAN 


The  VTOL  pitch  fan  is  utilized  in  the  XV5A  aircraft  and  was  tested  as 
mounted  in  the  aircraft  and  in  a  static  test  facility  (Figures  1  and 
41). 

fyj,  “  36  inches 


W  =  144  lb/sec  @  4074  RPM 

RFM  =  4074  maximum 

Nr  =  36 

Ns  =  90 

P  =  1.08 

r 
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Figure  41.  YTOL  Static  Test  Facility 


CRUISE  FAN 


The  cruise  fan  configuration  is  of  the  same  basic  design  as  the  VTOL  lift 
fan,  but  is  used  in  a  pod  type  installation  as  compared  to  a  wing  instal¬ 
lation. 


SINGLE-STAGE  SCALE  MODEL  COMPRESSOR  (Figures  42,  43,  44,  and  45 ) 

The  single-stage  scale  model  compressor  was  capable  of  having  a  variable 
rotor  -  stator  spacing  with  approximately  .5  inch  of  travel  in  the  stator 
row.  A  significant  characteristic  is  the  high  hub-tip  ratio  or  small 
ann’.lus  height,  which  is  in  keeping  with  the  basic  assumption  discussed 
in  Jection  One. 

=  26  inches 

=  .846 

=  11  lb/sec 
=  1780 
=  56 
=  48 
=  1.008 

DEVELOPMENT  VEHICLE 


The  data  presented  in  Figure  13  and  in  the  following  tables  under  develop¬ 
ment  vehicle  are  the  results  of  a  scale  model  fan  test  utilizing  a  fan 
configuration  of  proprietary  nature.  The  noise  data  is  not  classified  or 
proprietary. 
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Figure  42.  Single-Stage  Scale  Model  Compressor. 
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R^tor-Stator 
Spacing;  (Design) 


Annulus 

Height 


Figure  44.  Test  Configuration  of 
(reference  Table  1). 
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Bellmouth 


Model  Compressor 
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TEST  PROCEDURE  AND  EQUIPMENT 


FREE  FIELD  MEASUREMENT 


Free  field,  noise  measurement  of  the  various  scale  model  fan  and  compressor 
vehicles  were  made  in  generally  the  same  manner.  The  test  vehicles  were 
all  driven  by  a  25-horsepower  constant- speed  electric  motor  that  was  fed 
through  a  bank  of  reheostats  used  to  bring  the  motor  up  to  speed.  When 
steady  speed  was  attained,  a  microphone  was  placed  at  the  desired  location, 
and  the  sound  was  taped  for  a  minimum  of  15  seconds;  then  the  microphone 
was  moved  to  the  next  location.  This  procedure  was  the  same  for 
circumferential  surveys,  Figure  46,  arid  for  radial  traverses  of  the 
vehicle  annulus  over  the  rotor  face.  A  l/2-inch  microphone  was  used  with 
a  nose  cone  as  required.  The  various  systems  utilized  are  shown  in 
Figures  47  and  48. 

The  cruise  fan  and  XV5A  aircraft  free  field  measurements  were  made  with 
the  recording  system  shown  in  Figure  48.  The  fan  or  aircraft  was  brought 
to  the  desired  speed;  then  the  data  was  recorded  at  the  various  circum¬ 
ferential  locations,  Figures  49  and  50,  with  a  sufficient  recording  time 
for  data  analysis. 


VTOL  STATIC  TEST  FACILITY  MEASUREMENTS 


Measurements  made  of  the  VTOL  lift  and  pitch  fans  on  the  VTOL  static  test 
facility.  Figures  51  and  52,  were  done  using  a  multiple  channel  recorder 
with  17  microphones  positioned  around  the  fans.  The  applicable  operating 
conditions  were  provided  by  the  test  facility  operator  when  the  fans  were 
at  the  desired  RFM.  At  the  desired  operating  condition,  the  microphone 
outputs  were  recorded  on  the  multiple  channel  recorder  in  sets  of  three. 
The  instrumentation  is  shown  schematically  in  Figure  48. 


REVERBERATION  ROOM  MEASUREMENTS 


Reverberation  room  tests  of  the  treated  and  untreated  exhaust  louvers  were 
made  using  a  noise  source  shown  in  Figure  53 5  the  speaker  was  mounted 
in  the  plenum  chamber  to  which  the  louvers  were  attached.  The  vacuum  tube 
voltmeter  was  used  to  monitor  the  input  to  the  speaker,  and  the  oscillo¬ 
scope  was  used  to  insure  that  the  signal  was  not  distorted.  The  signal 
was  measured  using  a  one-third  octave  analyzer  and  level  recorder  for 
comparison  with  and  without  the  treated  louvers. 
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Figure  46.  Test  Layout  and  Microphone  Station  for  Noise 
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Figure  47  .  Data  Acquisition  System  -  Field  Test. 
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Figure  48  .  Data  Acquisition  System  -  Field  Test. 
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Figure  49.  Test  Layout  of  XV-5A  Research  Aircraft 
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Figure  52.  Test  Layout  for  Lift  Fan  Noise  Evaluation. 


81 


DATA  ANALYSIS 


DATA  REDUCTION 


Data  reduction  of  the  noise  tapes  made  on  the  various  vehicles  was  done 
using  the  two  systems  shown  schematically  in  Figure  5^*  The  spectrum 
analyzer  provided  an  analysis  at  a  10-cycle  bandwidth,  while  the  narrow 
band  analyzer  was  a  constant  6  percent  oandwidth.  The  6  percent  band¬ 
width  was  sufficient  for  reduction  of  the  first,  second,  and  occasionally 
the  third  harmonic  noise  level  of  all  the  vehicles  tested.  Noise  levels 
of  higher  harmonics  were  obtained  using  the  10-cycle  bandwidth  analyzer. 
The  6  percent  filter  was  not  used  when  the  harmonic  peak  was  not  at  least 
5  db  above  the  background  or  white  noise. 


SOUND  POWER  CALCULATIONS 


The  sound  power  levels  were  calculated  in  two  manners  corresponding  to  the 
circumferential  surveys  at  constant  radius  and  the  radial  measurements  in 
the  vehicle  annulus  near  the  rotor  face. 

The  circumferential  survey  sound  pressure  levels  were  put  into  a  computer 
program  designed  to  calculate  the  sound  power  level  based  on  three  fourths 
of  the  surface  area  of  a  sphere  having  a  radius  equal  to  the  measuring 
radius.  The  three-fourths  surface  area  is  a  compromise  between  total 
ground  reflection,  50  percent  surface  area,  and  total  ground  absorption, 
100  percent  surface  area. 

The  sound  power  level  was  calculated  from  the  radial  sound  pressure  levels 
in  the  vehicle  annulus  using  a  technique  developed  in  previous  compressor 
measurements  (reference  l4,  page  3*)  The  power  level  is  calculated  from 
the  following  equation: 


FWL 


where  A  is  the  annular  area  corresponding  to  each  radial  SPL  measurement. 
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Figure  54  .  Data  Reduction  Instrumentation. 
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TABLE  2.  TEST  DATA,  INLET  GUIDE  VANE  -  ROTOR  FAN 


|wf*fin  — nt  Location:  see  Figure  46  Fan  Speed:  1780  RMC 

Hf«eure»snt  fodlust  20  ft  _ Fan  Fundamental  Frequenci:  1395  c/a 


Sound  Preesure  Level 
(db  re  .0002  microbar) 

Sound  Pover^Level 
(db  re  10J  umtts) 

<>v>*n8leJ  0 

Harmonic''>s. 

D 

20 

50 

80 

no 

130 

150 

160 

121 

U3 

no 

108 

104 

9  In.  Inlet  duct.  0*  10V 
setting  and  exhaust  louver 
setting. 

Fundaaental 

Second 

Third 

Fourth 

Fifth 

70 

65 

73 

68 

62 

91 

80 

76 

70 

70 

86 

78 

78 

71 

68 

82 

71 

75 

61 

64 

83 

75 

72 

66 

68 

89 

72 

70 

79 

71 

87 

85 

77 

72 

73 

85 

78 

71 

7k 

66 

''^Angle,  0 

30 

50 

80 

no 

130 

150 

170 

180 

Harmon 1 c 

fundamental 

74 

84 

85 

80 

79 

85 

73 

81 

89 

U9 

Second 

66 

78 

72 

68 

76 

78 

70 

78 

84 

111 

9  In.  Inlet  duct.  0*  lev 

Third 

7* 

72 

78 

74 

71 

76 

68 

75 

73 

no 

setting  and  20*  exhaust 

Fourth 

73 

70 

73 

66 

68 

71 

65 

70 

78 

106 

louver  setting. 

Fifth 

63 

73 

64 

64 

71 

72 

63 

66 

76 

105 

^^Angle,  0 

■ 

30 

50 

80 

no 

130 

160 

180 

Harmonic''-^ 

Fundamental 

72 

82 

87 

79 

79 

76 

81 

85 

118 

Second 

67 

73 

75 

67 

82 

76 

79 

80 

112 

9  In.  Inlet  duct.  0*  1GV 

Third 

71 

68 

75 

75 

72 

71 

74 

80 

109 

setting  and  30*  exhaust 

Fourth 

68 

67 

72 

66 

68 

73 

72 

80 

106 

louver  setting. 

Fifth 

65 

73 

67 

66 

69 

72 

73 

78 

106 

''"\Angle,  0 

■I 

B 

20 

40 

60 

80 

100 

120 

140 

160 

180 

Harmonic'"'-^ 

Fundamental 

76 

91 

92 

86 

79 

83 

89 

89 

89 

89 

123 

Second 

71 

85 

76 

78 

71 

78 

71 

78 

75 

82 

ni 

9  In.  inlet  duct.  0*  IGV 

Third 

75 

78 

74 

69 

66 

73 

74 

74 

73 

76 

108 

setting  and  40*  exhaust 

Fourth 

68 

69 

70 

65 

65 

74 

71 

71 

74 

74 

104 

louver  setting. 

Fifth 

66 

70 

73 

68 

65 

72 

72 

73 

70 

73 

105 

'-'^Angle,  0 

0 

20 

40 

60 

80 

100 

120 

140 

160 

180 

Fundamental 

93 

104 

104 

95 

97 

98 

103 

99 

105 

98 

136 

Second 

100 

99 

100 

96 

90 

93 

89 

94 

96 

90 

131 

4  In.  Inlet  duct.  -27*  IGV 

Third 

84 

87 

91 

86 

86 

86 

87 

89 

87 

85 

123 

setting  and  0*  exhaust 

Fourth 

81 

82 

79 

76 

75 

74 

81 

85 

81 

77 

U5 

louver  setting. 

Fifth 

76 

79 

78 

74 

70 

76 

80 

80 

82 

74 

U3 

'^Angle,  9 

m 

B 

20 

40 

6o 

80 

100 

120 

l4o 

160 

180 

Harmon  1  c'''-^ 

Bfl 

Fundaaental 

87 

97 

98 

87 

82 

90 

96 

90 

94 

90 

129 

Second 

83 

63 

92 

83 

79 

81 

88 

89 

90 

80 

122 

9  In.  Inlet  duct.  -28*  IGV 

Third 

82 

79 

80 

83 

80 

?6 

82 

81 

81 

73 

U6 

setting  and  0*  exhaust 

Fourth 

76 

74 

77 

72 

73 

79 

76 

77 

74 

69 

in 

louver  setting. 

Fifth 

74 

76 

74 

68 

66 

66 

76 

76 

76 

68 

109 

'''-vAngle,  0 

mm 

B 

20 

40 

60 

80 

100 

120 

140 

160 

180 

Harmrnlc 

Fundamental 

80 

96 

91 

85 

84 

80 

93 

86 

84 

76 

125 

Second 

8o 

83 

85 

87 

81 

76 

82 

83 

80 

78 

U9 

9  in.  Inlet  duct,  -lb*  IGV 

Third 

79 

82 

83 

78 

72 

74 

82 

82 

80 

68 

U5 

setting  and  0*  exhaust 

Fourth 

68 

73 

78 

70 

63 

66 

76 

79 

74 

67 

no 

louver  setting. 

Fifth 

68 

-Jt- 

-22- 

68 

6? 

P 

Jj- 

67 

106 

85 


TABLE  3.  TEST  DATA,  ROTOR-STATOR  FAN 


Heaaureoent  Location:  Kf  Figure 
Hcuurcaen'  Radius:  10  ft 


r»n  Speed:  I78O  FfH 

ru  FundMhhthl  Frequency:  1070  c/» 


0  10  20  50  90  50  60  ro  SO  90  100  HO  120  150  190  150  160  176  100 


rundventel  81  79  88  92  85  90  8/  80  85  85  09  87  69  91  77  89  8r  90  87 

Second  75  76  77  8j  79  62  80  70  7  5  79  81  7i  79  82  81  76  80  77  76 

Third  68  72  77  79  81  78  80  78  77  70  77  7  2  77  ‘,5  80  72  77  -  77 

rourth  71  70  66  71  70  75  72  66  6}  66  6‘,  75  7*  72  75  76  65  - 

rifth  65  67  76  70  70  65  65  66  65  59  66  69  71  67  75  72  68  - 


rundearnua  8r  85  95  9»  86  92  95  86  82  86  90  90 

Second  82  82  85  81  88  8  5  85  78  78  80  85  77 

Third  75  81  8*  85  89  89  82  87  89  85  90  8o 

rourlh  75  71  68  71  68  75  70  69  72  72  80  7“ 

rifth  85  75  80  75  68  76  75  65  62  7o  60  78 


10  20  50  90  50  60  70  80  90  100  110  120  150  190  150  160 


Fund— ental 

Second 

Third 

Fourth 

Fifth 


75  00  82  80  84  86  84  f*7 


100 

110 

IPO 

1)0 

140 

150 

16© 

170 

90 

90 

/y 

91 

86 

90 

66 

92 

120 

55 

77 

da 

82 

59 

73 

00 

77 

11) 

90 

e? 

•74 

8^ 

80 

77 

84 

74 

116 

ao 

77 

TJ 

7) 

75 

70 

74 

104 

6 6 

75 

75 

ei 

61 

75 

75 

72 

105 

Fund— entsl 

8? 

76 

66 

92 

88 

82 

82 

05 

84 

70 

85 

85 

77 

80 

84 

79 

85 

Second 

80 

72 

77 

70 

80 

78 

75 

78 

78 

75 

69 

64 

64 

70 

75 

75 

66 

Third 

75 

72 

75 

77 

77 

72 

78 

72 

71 

70 

69 

59 

77 

75 

75 

60 

68 

Fourth 

70 

66 

75 

70 

75 

69 

68 

60 

65 

60 

58 

54 

61 

60 

56 

59 

56 

Fifth 

74 

78 

85 

77 

75 

76 

.5 

64 

65 

57 

58 

68 

65 

62 

65 

64 

62 

Fan  equipped  with  4  it;. 
Inlet  duet  sod  stator 
row.  No  exhaust  louvers. 


Angle,  6 

0 

10 

20 

50 

40 

50 

60 

70 

80 

90 

100 

110 

120 

1)0 

140 

150 

— 

Ha™onlc>sss^ 

Fund— entaj 

88 

80 

88 

91 

85 

88 

89 

82 

84 

86 

89 

68 

92 

66 

86 

90 

117 

Fan  equipped  with  4  In. 

Second 

85 

85 

89 

80 

85 

80 

78 

80 

79 

78 

79 

84 

84 

82 

89 

85 

112 

Inlet  duct  and  stator 

third 

76 

77 

75 

80 

75 

74 

80 

75 

74 

71 

80 

80 

76 

81 

60 

77 

107 

row.  Stator  row  spaced  1.0 

Fourth 

67 

72 

68 

79 

72 

74 

71 

72 

71 

75 

74 

n 

75 

75 

74 

79 

102 

In.  fro*  rotor  u  cc*  pared  to 

IHfiH 

74 

69 

75 

67 

71 

71 

b8 

64 

65 

67 

06 

66 

7- 

77 

74 

72 

100 

nonxal  spacing  of  .575  In. 

Fan  equipped  with  4  In. 
Inlet  duct,  stator  row,  and 
12  In.  acoustically  treated 
exhaust  louvers  (0*  louver 
a.*le). 


100 

no 

120 

1)0 

140 

150 

160 

170 

85 

81 

88 

92 

90 

91 

65 

99 

120 

62 

87 

66 

66 

87 

00 

00 

91 

llu 

Fan  equipped  with  4  lr.. 

88 

89 

66 

05 

84 

75 

79 

68 

U5 

Inlet  duct,  stator  row,  and 

71 

72 

75 

77 

75 

70 

79 

8)  1 

104 

standard  exhaust  louvers 

V> 

69 

75 

81 

79 

74 

75 

79  | 

106 

(0*  louver  angle) 

TABLE  4.  TEST  DATA,  ROTOR-STATOR  FAN 


Measurement  Location:  see  Figure  46 
Measurement  Radius:  20  ft 

Fan  Speed:  1780  RPM 

Fan  Fundamental  Frequency:  1070  c/s 

SSI  | 

Configuration 
(reference  Figure  39) 

'\Argie,  fl 

HanBonicv-u. 

B 

20 

40 

80 

100 

120 

140 

160 

117 

109 

112 

106 

106 

4  in.  inlet  duct.  Ho 
exit  duct  or  plug. 

Fundamental 

Second 

Third 

Fourth 

Fifth 

88 

72 

7** 

74 

84 

7e 

74 

72 

72 

70 

75 

72 

72 

82 

76 

78 

69 

66 

76 

72 

80 

68 

69 

80 

70 

73 

74 

73 

85 

76 

82 

72 

73 

84 

70 

75 

72 

75 

•^Angle,  0 

C 

20 

4o 

60 

80 

IOC 

120 

140 

160 

Harmonic^. 

Fundamental 

88 

81 

76 

81 

81 

83 

87 

90 

84 

120 

Second 

83 

74 

83 

77 

74 

74 

77 

82 

76 

114 

9  in.  inlet  duct.  No 

Third 

80 

78 

79 

80 

75 

77 

77 

84 

75 

114 

exit  duct  or  plug. 

Fourth 

80 

77 

74 

71 

70 

72 

75 

77 

72 

109 

Fifth 

78 

76 

72 

69 

66 

70 

72 

76 

74 

177 

v^Angle,  a 

■1 

20 

40 

60 

80 

100 

120 

140 

160 

Karmonlc^-^. 

Fundamental 

85 

82 

85 

80 

82 

80 

78 

80 

76 

116 

Second 

81 

83 

78 

71 

72 

71 

78 

82 

75 

113 

23-5  In.  inlet  duct. 

Third 

75 

76 

78 

76 

74 

80 

75 

62 

78 

113 

No  exit  duct  or  plug. 

Fourth 

70 

72 

72 

69 

66 

68 

71 

- 

72 

105 

Fifth 

84 

78 

70 

o7 

64 

67 

71 

75 

7C 

106 

'^sAngle, 0 

■I 

■ 

20 

4o 

60 

80 

100 

120 

140 

160 

H 

Harmonic''^ 

w‘ 

■ 

Fundamental 

81 

79 

74 

75 

76 

72 

82 

82 

78 

Second 

75 

76 

76 

74 

69 

72 

75 

76 

66 

9  in.  inlet  duct.  16  in. 

Third 

80 

78 

76 

73 

74 

73 

70 

76 

74 

■  ~ 

exit  duct  and  cylindrical 

Fourth 

78 

76 

74 

68 

64 

64 

75 

74 

71 

107 

exit  plug. 

Fifth 

77 

73 

72 

70 

64 

" 

69 

77 

72 

106 

■'^vAngle,  0 

mm 

20 

40 

60 

8o 

100 

120 

140 

160 

Harmon  lc''-^ 

Fundamental 

86 

80 

82 

79 

79 

68 

82 

78 

82 

115 

Second 

76 

75 

72 

74 

75 

69 

68 

77 

74 

109 

9  in.  inlet  duct.  16  in. 

Third 

82 

72 

77 

76 

70 

73 

71 

76 

74 

309 

exit  auct  and  conical 

Fourth 

80 

74 

73 

71 

67 

66 

74 

74 

69 

107 

exit  plug.  Fan  exit  duct 

Fifth 

79 

74 

72 

70 

63 

64 

70 

77 

74 

107 

area  with  plug  »  446  In. 2 

87 


TABLE  5.  TEST  DATA,  ROTOR-STATOR  FAN 


Sour-d  Pressure  level 
(db  re  .0002  iljcbir) 


0 

10 

20 

.30 

40 

50 

60 

TO 

so 

90  100 

110 

120 

130  140 

130  160  170 

180 

77 

80 

79 

78 

60 

75 

71 

73 

75 

73 

66 

T» 

(Z 

76 

71 

77 

73 

75 

116 

Fan  equipped  with 

64 

68 

68 

65 

70 

69 

06 

65 

60 

58 

62 

62 

(£ 

71 

65 

67 

60 

64 

65 

107 

4  in.  lr.let  duct. 

61 

59 

61 

64 

66 

6? 

66 

67 

61 

59 

56 

oO 

55 

5> 

60 

59 

64 

57 

62 

104 

Ho  stator  row  or 

56 

57 

62 

63 

56 

6^ 

65 

53 

58 

51 

56 

55 

62 

53 

55 

52 

52 

57 

57 

101 

exhaust  louvers. 

51 

57 

58 

57 

39 

57 

56 

51 

50 

51 

48 

48 

53 

58 

59 

51 

48 

50 

52 

96 

0 

10 

20 

30 

U) 

50 

60 

70 

so 

90  100 

no 

120  1^0 

1-0 

ro 

160  170 

160 

84 

82 

84 

60 

78 

73 

77 

77 

77 

75 

71 

72 

73 

80 

82 

72 

87 

82 

71 

12C 

Fan  equipped  with  4 

71 

75 

76 

76 

69 

67 

72 

69 

68 

67 

66 

ri 

72 

69 

75 

66 

75 

75 

67 

115 

m.  Inlet  duct  and 

66 

69 

69 

65 

66 

75 

69 

73 

72 

70 

7** 

71 

'’O 

t-6 

ofc 

68 

66 

0«* 

60 

112 

stater  row.  2  tat  or 

5e 

65 

64 

65 

66 

58 

61 

56 

55 

53 

57 

56 

66 

63 

64 

t>7 

63 

60 

56 

104 

row  spaced  1  ir.. 

61 

58 

60 

63 

63 

58 

57 

55 

54 

51 

53 

53 

59 

ol 

56 

6t> 

66 

57 

54 

101 

fro*  rotor  as  coc- 
pared  to  ncraal 
spacing  of  -575  in. 

0 
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80 
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84 

77 

78 

79 

so 

76 

71 

72 

72 

74 

67 

66 

66 

73 

79 

77 

70 

116 

Far*  equipped  with  4 

64 

67 

67 

66 

71 

67 

67 

59 

64 

54 

59 

66 

62 

71 

59 

61 

60 

58 

107 

In.  Inlet  duct. 

61 

56 

65 

61 

67 

66 

6* 

66 

65 

62 

66 

59 

64 

58 

56 

50 

59 

51 

105 

stator  row,  and  12 

57 

61 

68 

64 

66 

59 

61 

54 

55 

51 

56 

56 

60 

63 

58 

58 

58 

51 

102 

In.  acoustically 

74 

fl 

72 

66 

69 

69 

61 

62 
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TABLE  7.  TEST  DATA,  PITCH  FAN 
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APPENDIX  I.  JET  ENGINE  COMPRESSOR  NOISE 


PURPOSE 


The  purpose  of  this  work  is  to  review  investigations  on  the  generation 
of  noise  in  the  compressor  assembly  of  jet  engines  and  its  subsequent 
propagation.  Most  of  the  reported  material  is  not  new  but  has  been 
assembled  from  worhs  that  are  given  in  the  bibliography  section.  It 
has  been  edited  in  a  manner,  considered  by  this  author,  as  systematic 
and  logical  in  order.  The  problem  of  the  noise  of  a  compressor  assembly 
of  a  jet  engine  is  reduceable  to  a  problem  of  generation  and  propagation 
of  sound  waves.  If  classical  assumptions  are  acceptable  physically, 
this  type  of  investigation  turns  into  an  eigenvalue  problem  with  the 
wave  equation  as  its  field  equation.  Thus  the  previous  worx  of  others 
is  not  essential. 


Throughout  Appendix  I,  the  following  symbols  will  be  used  in  place  of 
those  listed  on  page  xv . 


(x,  y,  z) 
(r,  9,  z) 

(£’*•£> 
t  ,T 

c 

w 

V 


Cartesian  coordinates 

cylindrical  coordinates 

nondimens ional  cylindrical  coordinates 

time,  nondimens ional  time 

velocity  of  sound 

argular  velocity  of  rotor 

stream  velocity 


u 

0 

V  ri 
0 

6  <  ) 

Jn>  Yn 
n,  m,  J,  k 

N  ,  N 
r  s 

h 


velocity  of  pressure  disturbance 

v  /c 
z 

radii,  outer  arid  inner  cylinder 

velocity  potential 

Dirac  delta  function 

Bessel  functions 

summation  indexes 

number  of  blades,  rotor,  stator 
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Interest  in  the  study  of  compressor  noise  of  jet  engines  centers  aWut 
toe  problem  of  the  generation  of  sound  in  the  compressor  and  its 
subsequent  propagation  to  an  observer  somewhere  in  space.  The  gener¬ 
ation  of  noise  is  closely  associated  with  the  pressure  fluctuations 
inflicted  upon  the  medium  under  consideration  by  the  compressor  assembly. 
Hie  propagation  of  noise  is  thus  identical  to  the  propagation  of  the 
pressure  fluctuations.  If  one  is  interested  in  investigating  these 
problems  analytically,  he  would  have  to  describe  the  nature  of  prop¬ 
agation  of  the  pressure  variations  by  means  of  equations  called  field 
equations  representing  faithfully  this  phenomenon. 

The  formal  way  of  obtaining  the  field  equations  is  to  start  with  the 
fundamental  lavs  of  mechanics:  namely  (l)  the  conservation  of  matter, 

(2)  the  conservation  of  momentum,  and  (j)  the  conservation  of  energy. 
These  are  supplemented  by  the  thermodynamic  equations  of  state  for  the 
medium  under  consideration.  If,  however,  it  is  assumed  that  the  problem 
under  consideration  adheres  adequately  to  the  assumptions  of  classical 
acoustics,  the  above  yield  the  well-known  wave  equation 

-  h  4lS  -°  W 

as  the  field  equation  for  the  problem. 

To  complete  the  formulation  of  the  problem,  it  remains  to  impose 
boundary  and  initial  conditions  upon  it.  These  are  usually  dictated  by 
the  geometric  and  physical  aspects  of  each  individual  case.  For  a  jet 
engine,  for  example,  the  boundary  Is  the  core  of  the  engine  on  the  in¬ 
side  and  the  casing  on  the  outside.  On  the  boundary  it  is  reasonable  to 
assume  that  the  velocity  normal  to  the  surface  of  the  boundary  is  the 
same  as  that  of  the  boundary  itself.  Thus,  if  the  core  and  the  casing 
are  considered  rigid,  the  velocity  of  the  medium  in  the  normal  - to - the - 
surface  direction  is  zero. 

Another  condition  can  be  extracted  from  the  continuity  of  flow.  Since 
in  classical  acoustics  changes  of  the  state  of  medium  are  neglected, 
then  the  velocity  of  the  fluid  in  the  axial  direction  should  have  no 
jumps  or  discontinuities.  Therefore,  it  is  required  that  the  axial 
velocity  be  continuous.  The  last  condition  is  formulated  from  the 
action  of  the  blades  of  the  rotor  and  of  the  stator  on  the  flowing 
medium.  If  the  blades  of  the  rotor,  say,  are  equally  spaced  and  if  the 
rotor  is  rotating  at  some  angular  velocity,  w  ,  then  the  tangential 
component  of  the  velocity  should  suffer  a  discontinuity  every  time  a 
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blftJ*  passes  by.  Thus  this  condition  can  be  expressed  as  a  finite  dis 
continuity  of  the  magnitude  of  the  tangential  velocity.  Now,  xhe 
analytical  protlem  of  noise  generation  and  propagation  is  completely 
defined  by 


- 


1 


^0 


-  0 


«1  (rl}  '  U1  (ro}  m  0 

u  J  =  continuous  (2) 

s  =  0 

=  f  (r)  §  (  0  -  vt  ) 

B  =  0 


where  u  =  (up,  U2,  u^)  -■  grad  0  ,  (f  is  the  acoustic  potential,  and £\ 
is  a  Dirac  delta  function  indicating  a  Jump. 


The  solution  of  this  mathematical  problem  will  predict  how  the  pressure 
fluctuations  created  by  the  compressor  assembly  propagate  through  the 
Jet  engine  toward  the  intake  and  exhaust  of  it.  Once  the  pressure 
waves  find  themselves  at  the  open  ends  cf  the  jet  engines,  they  will 
disperse  through  the  surrounding  space.  This  is  called  radiation  of  the 
pressure  waves,  and  the  manner  by  which  they  radiate  will  determine  the 
intensity,  directivity  and  frequency  by  which  they  reach  an  observer. 
Here,  again,  some  clarifications  of  the  boundary  conditions  of  the 
radiating  surface  are  necessary.  One  simple  way  of  looking  at  this  part 
of  the  problem  is  to  imagine  that  the  whole  space  is  separated  into  two 
parts  by  a  surface  containing  the  radiating  faces.  The  radiating  faces 
(in  this  case  the  open  end  of  the  engine)  are  characterized  by  the 
pressure  fluctuating  as  predicted  by  the  propagation  part  of  the  solu¬ 
tion;  the  rest  of  the  surface  is  rigid.  This  is  all  that  is  necessary 
to  calculate  the  radiation  field  completely,  outside  of  a  description  of 
the  medium  itself. 


It  should  be  pointed  out  that,  while  the  above  statements  contain  all 
the  fundamental  parts  necessary  for  the  study  cf  noise  generation  and 
propagation  by  the  compressor  assembly  of  a  jet  engine,  it  is  also  a 
simplification  of  the  overall  problem,  as  it  will  be  shown  subse¬ 
quently. 

A  schematic  representation  of  the  problem  of  generation,  propagation 
and  radiation  of  noise  is  given  in  Figure  55» 
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Figure  55  •  Evolution  of  the  Noise  Problem. 


With  this  brief  introduction  to  the  fundamental  problem,  one  now  can 
undertake  the  formulation  and  solution  of  the  jet  engine  noise  in  a 
systematic  manner.  The  basic  approach  is  rather  fixed,  as  indicated  in 
Figure  55.  There  is,  however,  quite  a  choice  in  the  order.  Sane 
authors  prefer  to  start  with  the  effect  of  the  rotor  and  of  the  rotor  - 
stator  Interaction  on  the  solution  of  the  field  equations.  Others 
choose  to  formulate  the  eigenvalue  problem  first  and  obtain  its 
solution,  then  inject  the  rotor  and  rotor  -  stator  effects  to  narrow  the 
solution  to  the  specific  problem  at  hand.  The  final  results  should  be 
the  same  regardless  of  the  choice.  In  this  review  the  following  subjects 
shall  be  discussed. 

(1)  The  formulation  of  the  eigenvalue  problem  from  the  physical 
and  geometric  aspects  of  the  jet  engine. 

(2)  The  solution  and  its  meaning.  The  influence  of  the  rotor 
effects  and  rotor  -  stator  interaction  on  the  critical  speed. 

(5)  Important  aspects  of  the  obtained  solution.  A  discussion  of 

the  design  parameters  which  enter  into  the  significant  results 
of  the  solution  and  if  possible  their  effect  on  noise 
characteristics . 

The  radiation  of  the  acoustic  waves  into  surrounding  space  will  not  be 
discussed  in  thic  report  since  it  will  be  covered  extensively  in  sub¬ 
sequent  work. 


THE  MATHEMATICAL  PROBLEM 


The  first  commitment  made  here  is  adherence  to  the  classical  acoustics. 
The  meaning  of  this  assumption  is  that  the  behavior  of  the  pressure 
fluctuations  is  described  by  the  classical  wave  equation.  Inherently 
this  equation  of  acoustic  waves  implies  that  the  pressure  fluctuations 
are  very  small.  This  is  most  likely  an  incorrect  assumption,  but  the 
price  for  a  more  accurate  description  is  the  introduction  of  nonlinear 
terms  in  the  field  equations,  which  mean  an  extremely  difficult  problem 
from  the  very  start.  Thus,  it  is  probably  wiser  to  accept  the  lack  of 
proper  description  at  this  point  for  the  sake  of  accessibility  to  a 
formal  solution.  Therefore,  the  field  equation  employed  for  this  kind 
of  work  is 

V20-  -0 

c  dt 

where  0  is  the  acoustic  potential,  or  velocity  potential  as  sometimes 
called. 
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The  second  simplification  is  on  the  geometry  of  the  Jet  engine.  The 
core  and  the  enclosure  of  the  engine  are  approximated  by  tvo  concentric 
cylinders.  The  region  of  interest  to  this  problem  is  the  annulus 
contained  between  the  inner  and  outer  cylinders.  Furthermore,  it  is 
assumed  that  the  annulus  region  extends  from  -*»  to  +e o  in  the  axial 
direction.  The  exaggeration  is  imposed  in  order  to  make  the  subsequent 
mathematical  problem  more  accessible.  It  means  that  it  is  not  necessary 
to  have  any  boundary  conditions  at  some  finite  length  from  the 
compressor  assembly  to  worry  about.  The  introduction  of  this  simpli¬ 
fication  has  an  effect  on  the  magnitude  of  the  reflected  wave.  But  it 
has  been  pointed  out  that  this  influence  may  be  neglected  without  great 
sacrifice  to  the  practical  aspects  of  the  problem. 

Next  it  is  assumed  that  the  surfaces  (walls)  of  the  concentric 
cylinders  are  rigid.  This  implies-  that,  since  the  walls  are  nan- 
deformable,  the  mechanical  vibration  of  the  care  and  the  enclosure  of 
the  Jet  engine  are  neglected.  Experts  claim  that  the  coupling  between 
the  acoustic  pressure  waves  and  the  characteristics  of  the  surrounding 
structures  is  not  of  importance  in  the  problem  under  consideration. 

Figure  56  shows  what  the  Jet  engine  looks  like  after  the-  above  simpli¬ 
fications  have  been  Introduced. 


Now  introduce  a  cylindrical  coordinate  system  (r,  0,  a)  as  shown  in 
Figure  56.  Then  the  problem  put  in  mathematical  form  becomes: 
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0 


(4) 


U1  (r0>  ■  U1  (rl>  * 

|u2  (»  =  0)  |  *  f(r)  5  (e 

(a  =  0)  =  continuous 


-  vt) 


(5) 

(6) 


where  u^,  u^,  u^  are  the  velocity  components  in  the  r,  0  and  z  direc¬ 
tions  respectively.  Equations  (4)  through  (6)  are  the  boundary  condi¬ 
tions  of  the  problem.  Equation  (4)  is  a  consequence  of  the  rigidity  of 
the  walls.  Equation  (5)  expresses  the  fact  that  the  tangential  component 
of  the  velocity  varies  with  the  radius  and  that  its  magnitude  contains 
finite  discontinuities  due  to  the  chopping  of  the  blades.  Equation  (6) 
states  that  the  flov  in  the  axial  direction  is  continuous.  Condition 
(5)  further  states  dependence  of  the  velocity  on  the  angular  velocity 
of  the  blades  and  also  on  the  position  of  the  point.  Some  authors 
prefer  to  use  different  boundary  conditions  instead  of  (5) •  For 
example,  they  assume  that  the  pressure  distribution  at  the  plane  z  -  0 
is  known. 


p  (r,  0,  0,  t)  =  f  (r,  0,  0,  t) 


(7) 


This  is  not  unreasonable  since  one  can  measure  the  pressure  at  some 
convenient  plane  and  vise  it  as  a  datum.  This,  however,  commits  the 
problem  immediately  to  numerical  method*  and  also  divorces  the 
analytical  problem  completely  from  the  design  of  the  compressor.  The 
choice  here  is  mostly  dictated  by  the  experimental  approach  one  wishes 
to  use.  In  this  review  the  farmer  formulation  is  preferred  due  to  its 
possible  tie  with  compressor -blade  design. 

Introducing  the  nondimens i on al  coordinates  (  (^  ,  0,  £  )  such  that 


r  =  r 


0=0 


Z  =  r 


t  =  r. 


(8) 


Equation  (3)  expressed  in  the  cylindrical  coordinates  (K,  0,  f  )  and  T 
becomes  ^  ^ 

(9) 


1  0 _  ,1  0 

i  t  r  a 


0 


*■ 


ff# 
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Nov  assume  that  0  ®  0  (  £  ,  0,^  j  T  )  it'  of  the  farm 

0  -  R  (£)  *  (£)  0(6  -  "'O  (10) 

vhere  the (0  =  0)  (0  -  vt)  is  suggested  by  erudition  (5).  Further¬ 
more,  the  equal  spacing  of  all  blades  and  the  angilar  velocity  of  the 
rotor  suggest  that(0)  is  periodic,  and  of  the  font 


©■ 


la  (8  -  vt)  . 


therefore, 


4  -»<£j  *(^).ln  (e-wi  >  . 

Substituting  (ll)  into  (9),  one  obtains 

« <£>  -  v1^ *  ***£ 

K  ^3^  CnJn(  Xnm  T)  *  DnYn  (Anm£) 

2  2  2  2 
Xnm  'nw-» 

and  the  velocity  potential  has  the  farm 

0(£,  *,£/n  =  [v1^  +  V"^)  • 

+  DnYn(Xnm^)j  .  ^  -  w  T  )  j  • 


(11) 


(12) 


lCnJn 


a?) 


Nov  the  boundary  conditions  (4)  through  (6)  are  utilized  to  determine 
the  arbitrary  constants  introduced  into  the  analytical  problem.  The 
noadimens lonal  variables  introduced  by  equation  (8)  are  used.  Equation 
(M  yields 

0  (li#  0#£  f  T )  =  0  ;  0  (l,  0,^  ,  T )  =  0  ;  h  - 
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or 


Vi  (hXm>  +  Vi  (h\m>  *  0 
V  +  V  <Xm>  '  0 


Solving  for  Cq  and  Dq,  one  obtains  the  characteristic  equation 

J'  (\  )  Y'  (hX  )  -  J'  (h\  )  Y'  (X  )  =  0  ,  (14) 

n  ,vnm  n  Vim7  n  /xnin  n  nm  * 


whose  solution  yields  the  values  of  X  nm  which  satisfy  boundary  condi¬ 
tions  (4).  Note  again  that 


\  2  2  2  2 
An  m  =  n  v  -  m 


or  m 


2  2  \2 
lv  -  A  mn  . 


(15) 


Equation  (6)  requires  U_  ■  continuous  or 

J  «  -  0 


f  If..  . 


Of 

Substituting  (lj)  into  (l6),one  finds  that 


A  =  -B 
m  m 


Condition  (5)  can  be  approached  in  the  following  manner: 
vortex  of  intensity  p(r)  created  by  a  blade;  then 


(16) 

(17) 

consider  a 


Hr)  =  j)  u2ds  -  ro  J  |  u2  £ 


dd 


By  equation  (5)  this  becomes 

r<£>  -  '<£>. 

Thus  equation  (5)  can  be  written  as 

(  £  »e»v,  r  ) 


.  -Oil 


5  (  e  -  vr  )  ; 


then 


and  since  there  are  N  blades, 

u2  (£>3  *  0  *  r  ) 


5 

=  2 

j=i 


^"j(  isl  u(o  +  —  n  J 

rQ£  UV  N 


-  w  T  ) 


Introducing  the  potential  function  0,  one  obtains 

N 


|J_ 

3  e 


.0  *  r(£>  Z  +  -  wr> 


(18) 


under  the  assumption  of  identical  blades,  thus  equal  vortex  intensity 
p  ^  y  ^  Now,  vising  equations  (l3)  and  (17), 

b  *  (19) 

|||  .IE21"  Anmela(0-ur)  |cnJn(Xcm^)  +  DnYn(Xnm^jl. 


Then  by  Fourier -Bessel  expansion,  (l8)  and  (19)  yield 


nm 


-  rsnr-J^  £r<£>  *  <£>  ^ 


(20) 


where 


fi<£>  -  CA<Xmn£>  +  D»fn(k“£> 

is  the  normalized  form. 

Thus  the  completely  defined  solution,  <Jf  for  this  problem  is 

0  (£,#.£,r).  1  1  A^111  {(8-’'r  >  ♦AnV-X2^^ 

n=KN  m=l  v 


(21) 


where  K  =  ^  =  integer.  This  obviously  is  a  restriction  on  n,  which 
comes  from  mathematical  considerations;  but,  as  will  be  shown  later,  it 
is  physically  necessary  also.  It  merely  lays  that  the  index  n  has  to 
be  a  multiple  of  the  number  of  blades,  N. 


no 


It  is  possible  now  to  draw  some  conclusions  without  considering  the 
rotor  -  stator  assembly  any  further.  Examining  equation  (21), it  is 
evident  that  the  critical  component  of  the  solution  is  the  exponential 
term, 

ein  0-wT)  +J n2w2  -X2nm 

and  in  particular  the  coefficient  of  %  .  If  the  argument  of  the 
radical  is  a  positive  number,  then  the  solution  in  z-direction  will  be 
sinusoidal  in  nature;  thus  the  0  solution,  and  therefore  the  acoustic 
waves,  will  be  propagated  without  attenuation  in  z.  If  the  argument  of 


the  radical  is  negative, 
thus  the  acoustic  waves 

then  the  z-solution  is  exponentially  decaying; 
are  attenuating.  Therefore,  if 

n 

propagation  in  a 

(22) 

^  X  nm 

w< - 

n 

decay  in  a  1 

(23) 

and 

n 

w  * 

critical  speed. 

(24) 

This  w  0  is  quite  often  referred  to  as  the  cutoff  speed.  It  depends, 
so  far,  on  the  eigenvalues  of  the  wave  equation,  which  were  imposed  by 
the  rigidity  of  the  boundary.  It  will  be  shown,  in  the  next  paragraphs, 
that  the  critical  speed  depends  on  the  interaction  (or  interference)  of 
the  rotor  -  stator  system,  the  tip  Mach  number,  the  spin  velocity,  etc. 


THE  CRITICAL  SPEED 


In  the  previous  section  the  motion  of  acoustic  waves  in  the  jet  engine 
has  been  discussed  and  its  dependence  on  such  variables  as  the  geometric 
shape  of  the  engine  ( through  X  nm) .  the  speed  of  the  rotor  w,  and  the 
number  of  blades  (through  n  =  KN)  has  been  indicated.  The  influence  of 
the  rotor  -  stator  assembly  ^nd.  of  the  flow  velocity  on  the  critical 
parameters  (Xma,  w,  n,  F  ( L  )j  and  subsequently  on  the  critical  speed 
v  °  will  be  undertaken  in  €1118  section. 


Ill 


THE  ROTOR  EFFECTS 

The  formulation  of  the  problem  in  the  previous  section  it  essentially 
the  problem  of  generation  of  acoustic  waves  by  a  rotor.  It  remains  now 
to  justify  some  of  the  assumptions  made  there.  At  first  the  pressure 
fluctuations  are  created  by  the  rotating  blades.  Therefore!,  the  pres¬ 
sure  field  is  spinning  around  with  the  same  angular  velocity  as  the 
rotor.  Hence  if  one  expressed  the  pressure  with  respect  to  a  blade, 
the  pressure  would  be  expressed  in  terms  of  its  position  wilh  respect 
to  that  blade.  If,  however,  the  reference  frame  is  stationary,  then 
the  0  component  of  the  pressure  will  have  to  be  corrected  by  wt,  the 
arc  traveled  by  the  blade  in  time  t.  Thus  the  dependence  of  the  pres¬ 
sure  field  on  0  and  t  has  the  general  form 

P  (  e,  t)  =  P  (  0  -  wt). 

This  is  consistent  with  the  form  of  traveling  waves. 

If  the  blades  are  assumed  to  be  equally  spaced  and  identical  Ln  design, 
then  the  pressure  fluctuations  will  be  periodic  in  nature  with  a  period 
1  W  being  the  number  of  blades.  Thus  the  pressure  functioi  has  the 
Nv 
form 

P(  0,  t)  =  P  (n  (  e-wt)J  , 

which  is  the  same  as  obtained  previously,  where  the  index  oi  summation 
is  n  s  kN,  k  =  integer,  rather  than  n  =  1,  2,  3  ....  As  foi  the  crit¬ 
ical  speed  w*  ,  the  previous  remarks  remain  unchanged,  for  the  "rotor 
only"  case. 


THE  ROTOR  -  STATOR  INTERFERENCE 


The  interaction  of  the  rotor  -  stator  assembly  is  caused  by  the  cutting 
of  rotor  waves  by  the  stator  blades,  the  cutting  of  the  stater  wakes  by 
the  rotor,  and  interference  on  rotating  pressure  patterns.  'Haile  their 
relative  influence  depends  on  the  geometry  of  the  blades,  thrix*  rela¬ 
tive  distance,  the  aerodynamic  aspects  of  flow,  and  whether  one  deals 
with  viscous  wakes  or  nonviscous  vorticity  sheds,  they  a’l  have  one 
common  property:  the  same  periodic  recurrence.  Therefore,  the 
periodic  behavior  of  these  interactions  will  be  discussed  first. 

Consider  a  rotor  and  a  stator,  concentric  to  each  other,  the  rotor  having 

N  blades  and  the  stator  having  N  blades,  both  being  equally  spaced 
I*  s 
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over  a  2  *  Interval.  The  periodic  nature  of  the  pressure  field  in  6 
and  t  is  expressed  by  the  relationship 

p  (e.t)  »  p  (e  +5,  t  +  l)  (25) 

The  intervals  6  and  t  can  be  obtained  as  follows:  Let,  at  time  t  ,  a 
blade  of  the  rotor,  say,  number  1,  be  located  opposite  a  blade  of  the 
stator,  say,  also  number  1.  If  Nr  t  Ns,  then  the  rest  of  the  blades  are 
not  centered  with  respect  to  their  rotor  or  stator  counterparts.  Since 
the  rotor  is  rotating  with  angular  velocity  w,  it  will  take 

t  =  ,*  A  0  =  -  §-£  f or  N  >N  (26) 

v  N  N  s  v 

s  r 


time  to  have  the  next  pair  of  rotor  and  stator  blades  centered  with 
respect  to  each  other.  This  relatively  equivalent  position  of  the 
stator  and  rotor  position  is  repeated  at  intervals  of  time  of  intergal 
multiples  of  A  9  .  Thus 

w 


(2?) 


In  the  period  of  time  t,  the  pressure  pattern  has  traveled  along  from 
stator  blade  1  to  a  neighboring  stator  blade,  say,  2.  The  arc  interval 
is  obviously  2  it  .  As  the  rotor  rotates  by  another  A  9,  the  pressure 
N 

s 

pattern  shifts  by  another  2  rt  .  Therefore, 


9  = 


2  n 


(28) 


The  angular  velocity  with  which  the  pressure  field  swirls  around  is  now 


N  N 

r  _  _r_ 

N  -N  W  "  AN 

r  s 


(29) 


2  it 

Next  it  is  realized  that  the  period  on  9  is  ^  ;  thus,  in  view  of  the 
comment  under  "The  Rotor  Effects",  one  can  express  the  pressure  field  as 


P  (e,t)  =\jj(e-M  t)  f  (t)  (30) 

where  f  (t)  is  a  time  function  due  to  the  fact  that  6  and  t  do  not 
rotate.  However,  since  the  frequency  on  time  is  w  ,  f  (t) 

9 
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may  be  expanded  in  series 


(t)  .  f  e  *2i«J  f  * 


(3D 


and  consequently  (30),  after  similar  expansion,  yields 


P  (6,t)  = 


£  ikM(e-vt)  -2jtiJ-r-  t 

J,k  6  6  0 


(32) 


or 


P  (0,t)  .  X  e'15  <“■  +  ¥  J)  ‘ 


(33) 


Examination  of  equations 
(l)  The  pressure 


(32)  and  (33)  yields  the  following  information: 
field  rotates  with  angular  velocity 


(34) 


(2) 


in  distinction  to  w  for  the  pure  rotor  case. 


The  general  expression  of  the  harmonics  generated  is  given 

by 


(kZsN 


2  « 
6 


•j)  =  w  kN  +  J 


N  N 


s 


(35) 


Equations  (32)  and  (33)  define  the  nature  of  the  acoustic  potential  in 
6  and  t  for  the  case  of  rotor  -  stator  interaction.  The  wave  propagation 
for  this  case  is  identical  to  the  problem  solved  previously  except  for 
the  (0)  component,  which  assumes  the  form  given  by  equations  (32)  and 
(33)*  Thus  the  potential  is  given  by 


(jt  -  r  c 
r  o 


Z  2  2 

k  j  m 


(36) 
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is  given  by 


where  row  the  coefficient  of 


+  J 


1/2 


with  the  restriction  on  n  ■  k  .  N,  n  =  integer. 


(??) 


The  critical  speed,  in  general,  will  be 

„  •  -  Xnm  A  H 

run  kNM  +  J  N  N 
r  °  r  s 


(58) 


Assuming  J  =  0,  since  it  appears  that  the  contribution  of  the  J's  will 
be  in  high  harmonics,  one  obtains 


w 

nm 


A  N  /  X  £aa\ 

N  '  n  ' 

r 


(59) 


Equation  (59)  compared  with  equation  (24)  indicates  that  the  critical 
speed  changes  by  a  factor  of  AN/Nj.,  with  frequencies  of  the  interaction 
between  rotor  and  stator  being  harmonics  of  Nr  w  (of  the  pure  rotor 
case),  for  the  case  of  J  =  0. 


THE  EFFECT  OF  CONSTANT  SUBSONIC  FLOW 


So  far,  the  problem  under  consideration  has  been  for  sound  propagation 
in  a  medium  standing  still.  The  case  of  the  stream  moving  with  a 
velocity  V  ,/  V  /  <C,  subsonic  and  constant  in  magnitude  will  be 
developed.  The  same  result  may  be  obtained  if  one  considers  the  medium 


at  rest  and  the  sound  source  moving 
equation  to  start  with  is  the  same. 


and  the  position  of  the  source  is  ( 
ducing  the  new  variables, 

x  = 

y  = 

5  = 

t  = 


with  ■ 

-  V  velocity. 

The 

field 

a2# 

dt2 

ii 

o 

>• 

i*  (t). 

,  y*  (t),  z* 

(t)  : 

) .  Intro- 

X  -  X 

*  (t) 

y  -  y 

*  (t) 

m 

z  -  z 

*  (t) 

t  . 


1 

i 

j 

i 
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The  transformation  introduces  the  change  in  0: 

0  (x,  y,  z,  t)  =  0^  ♦  x*  (t),  y  «■  y*  (t),  z  ♦  z*  (t );  1 j  .  (41) 

Substituting  into  the  wave  equation,  one  obtains  the  proper  form  of  the 
field  equation  for  a  moving  medium: 

V2  0  -  V  -  (iV)(v-V)  0  +  (  ^*7)4  =°- 

c  dt  c  L  * 

(42) 


Uniform  Axial  Flow 


This  case  is  characterized  by 


■[0 


,  o,  V 


zj 


v  =  constant . 
z 


Substituting  equation  (4})  into  (42)  yields 


2  . 

V  0 


.  k.  £& 

c2  St2 


2v, 

2 


x2  u  v  2 

o  0  _  z 

dtdz  c2 


d2  0 
TT2- 
o  z 


=  0 


(43) 


(44) 


Introduce  now  the  variables  in  cylindrical  coordinates 

C±T 


r  =  r 


C  ! 


z  =  r 


.2 


0=0 


t  =  r  X  ,  p.  i 

o  c  c 


(45) 


and  the  field  equation  assumes  the  form 

\2.  j.  „  „  ^2  0 


y20  .  +  -M 

V  i-  l 


(46) 
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The  solution  to  equation  (46)  is  obtained,  as  shown  previously,  by  the 
method  of  separation  of  variables,  and  is  of  the  same  form  as  equation 
(21).  The  difference  is  in  z,  where  the  factor  is  included 

t  ^  v/l-P 

in£  andT. 

Thus 

0  (£,  ,T)  .  eln  [(e  ‘  wT)  (Xnm£)  (47) 


and 


w  = 


Jl^2 


&  run  =  n<  - 


it 

17  w  V 

Xu  m  . 

£ 

wp 

l 

_U-  P2/ 

n2(l-p2)  . 

1-  e2J 

m 


The  critical  speed  is  obtained  by  setting  01 nm 
(26),  one  gets 


-  O 

W 

nm 


X 


nm 


or 


;  •  J 1-  p2  A- . 

nm  v  n 


=  .0j  thus, from  equation 

(49) 

(50) 


This  indicates  that  the  critical  speed  or  cutoff  speed  decreases  as  vz 
or  p  increases.  Similarly,  one  can  obtain  the  effect  of  moving  medium 
on  the  rotor  -  stator  system.  The  critical  or  cutoff  speed  in  the  case 
of  a  one-stage  compressor  becomes 


w  0 
nm 


\/l-  P 
kNr  +  j 


N  N 
s _ r 

A  N 


(51) 


in  contradistinction  to  equation  (58).  For  more  complicated  cases,  such 
as  the  inclusion  of  swirl  velocity,  one  proceeds  in  a  similar  way  by 
introducing  a  set  of  variables,  as  in  equation  (44),  where 
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e  7  T 


(52) 


@ 

where  v 


r 


o 


2 


is  the  swirl  velocity. 


7  = 


# 

w  r 

o 


c 
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APPENDIX  II.  THE  INCOMPRESSIBLE  PRESSURE  FIELD  OF  A  FLAT  PLATE  CASCADE 


Weinig  (reference  20)  concluded  that  the  incompressible,  inviscid  flow 
performance  of  any  cascade  could  be  described  in  terms  of  an  "equivalent 
flat  plate  cascade".  Weinig  and  many  others,  ineluding  the  present 
author  (references  12  and  13),  have  developed  techniques  for  quantitative 
demonstration  of  the  equivalence  between  real  cascades  and  flat  plate 
cascades.  With  the  help  of  modern  computing  machinery,  come  of  these 
methods  could  be  used  to  evaluate  pressure  fields  for  real  cascades.  For 
qualitative  evaluation,  however,  consideration  of  the  flat  plate  cascade 
can  be  useful. 

The  development  of  reference  20  shows  that  the  field  of  velocity  com¬ 
ponents  of  a  flat  plate  cascade  (Figure  57)  may  be  described  parametrically 
in  terms  of  a  "picture  plane",  Figure  58. 


Suitable  algebraic  expressions  are  available  relating  picture  plane 
parameters  h,(  ,  £  to  the  physical  plane  geometric  parameters  u  and 

B,  and  the  aefoayngmic  variables  £  and  The  identification  of 

values  of  the  parameter  corresponding  to  physical  plane  points  on  a 
desired  contour,  such  as  a  vertical  line  at  a  specified  distance  from 
the  cascade,  is  a  formidable  job  for  a  desk  calculation,  but  can  be  done 
expeditiously  with  the  help  of  the  IBM  709*+  computer.  Some  representative 
results  are  shown  on  Figure  2.  It  is  clear  that  the  magnitude  of  the 
velocity  fluctuations  seen  by  an  observer  moving  parallel  to  the  cascade 
axis  diminishes  very  rapidly  with  distance,  either  upstream  or  downstream. 

An  analytical  expression  for  the  decay  rate  may  be  deduced  from  equations 
(53)  and  (5*+)  .  In  the  region  far  upstream  frou  the  cascade,  equation 
(5*+)  is  effectively 


f  +  if  =  2T-  1"  (£♦  1  -  ih)  . 
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Figure  57.  Velocity  components  of  a  flat  plate  cascade 
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If  the  vector  distance  of  the  parameter  l %  from  the  image  of  upstream 
infinity,  (-1  +  ih),  is  expressed  in  polar  form,  r  (e)*®,  then 

In  (£+  1  -  ih)  =  In  r  +  i  9  (56) 

so  that  the  distance  between  two  points  (x^,  y^)  and  (x2,  y2)  is 


Then,  if  the  two  points  have  the  same  ordinate,  (y2  =  y^). 


Equation  (53)>  conveniently  expressed  in  logarithmic  form,  can  be 
developed  as  a  Taylor  series  based  on  the  Image  of  upstream  infinity: 

Wsln(ux-i  uy)=ln  c  +  In  (£-  £g)  ‘  1“  (£-£L) 

(£+  1  -  ih)  +  0  (£+  1  -  ih)2. 

(59) 

Neglecting  the  quadratic  and  higher  order  terms,  and  using  the  polar 
form  for  the  vector  distance,  as  in  (56), 

W  -W^  =  c'  r  e  19  .  (60) 

It  is  also  convenient  to  use  a  polar  form  for  the  velocity  vectors 

u  -  i  u  -  v  e 
x  y 

W  =  In  v  -  i  0 
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How,  applying  (60)  and  (6l)  to  a  point  far  upstream  from  the  cascade, 
where  the  local  velocity  magnitude  is  not  much  different  from  that  at 
infinity, 

In  i - i  (0  -  )  =  In  (1  +  -  ■  )  -  1  w  -  ) 


7  -  V 

OB 

~ - 1 


*  P„) 


=  r  Qte  (c 1 )  cos  0  -  Im  (c ' )  sin  9J 

+  in  £(Re  (c')  sin  0  +  Im  (c')  cos  sj - (62 ) 

Then  the  real  part  of  (62)  may  be  applied  to  the  two  points  (x. ,  y  )  and 
(x2,  y2)  8ame  ordinates  (y2=y1): 


v2  -  va 

vmt 


V1  ’  V« 


(63) 


Now,  substituting  (58)  into  (63), 
f2 


Y1  -  Vo* 


=  exp  I  2k 


-  V-  (x  -  X  ) 

2k  -= - =■ 

s 


If  (x2,  y2)  is  farther  from  the  cascade  than  (x1 ,  y1 ),  the  exponential 
is  negative,  demonstrating  a  rapid  decay  rate.  'This  rate  is  approxi¬ 
mately  a  factor  of  530  for  each  blade  spacing  away  from  the  cascade, 
which  Is  consistent  with  Figure  2.  The  rate  is  also  approximately 
55  db  for  each  blade  spacing,  and  appears  to  be  quite  close  to  Figure 
14  of  reference  18. 


The  analysis  downstream  of  the  cascade  is  similar,  except  that  the  log¬ 
arithm  in  (55)  is  negative,  which  leads  to  a  sign  change  in  the  exponent 
of  (64). 


123 


APPENDIX  III 


SUBSONIC  COMPRESSIBLE  FLOW  INFLUENCE  ON  THE 


PRESSURE  FIELD  OF  A  CASCADE 


In  the  stucy  of  noise  generation  and  propagation  from  the  compressors 
and  fans  of  aircraft  jet  engines,  one  of  the  mechanisms  investigated 
has  been  the  inviscid  pressure  field  associated  with  an  isolated  cascade. 
Reference  l8  •gfives  a  treatment  of  this  feature  which  makes  it  appear  to 
be  fundamentally  dependent  on  cylindrical  geometry.  Actually,  the  sub¬ 
sonic  results  appear  to  be  primarily  two-dimensional,  as  can  be  seen  in 
the  following  analysis  based  on  the  Karman-Tsien  approach. 

A  logarithmic  hodograph  plane,  Figure  59 >  is  the  basis  for  the  quasi- 
conformal  transformation  relating  the  complex  velocity  potential,  the 
complex  velocity,  and  the  physical  plane  coordinates,  Figure  60,  far 
away  from  a  cascade.  In  this  plane 

£=  ofe  -  i&o  (65) 

represents  the  flow  conditions  without  the  disturbances  due  to  indi¬ 
vidual  blades. 

£=  <u  -  <o„,  -  i  (  X  -  )  (66) 

' 

represents  the  disturbance  velocity,  which,  if  related  to  physical 
plane  coordinates,  produces  the  sound  sensation.  The  compressible 
velocity  function,  oj  ,  is  related  to  the  ratio  of  the  actual  velocity 
to  the  undisturbed  velocity,  X_, with  the  Mach  number  M  of  the  undis¬ 
turbed  flow  as  parameter,  by  ^ 

2  fer)  m~/xA  -  ^  (67) 

1  +  /l  +  M^2  /  Jl  -  M_S 

A  suitable  complex  velocity  potential  is 
oj  =  0  +  i  x 

=  T.  exp  (  ^  )  /  2  n  In  (£')  +  0  (  £').  (68) 

In  the  regime  of  interest,  the  terms  0  (  £')  may  be  neglected  in  compar¬ 
ison  with  the  logarithm. 
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Figure  59.  Logarithmic  liodograph  Plane 
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The  transformation  from  hodograph  plane  to  physical  plane  Is  defined 
by 

[,«,<£' )a£]  /(*£*>  (69) 

-  l/k  [b  exp  (£  )  exp  (2^o)  d  /  (2 *£). 

The  constant  B  Is  the  same  in  equations  (68)  and  (69)  and  is  to  be 
determined  so  that  one  circumnavigation  of  too  In  the  hodograph  plane 
corresponds  exactly  to  one  movement  across  a  blade  space  3  parallel  to 
the  cascade.  This  condition  is  satisfied  if  the  real  and  imaginary 
points  of  B  are  respectively 


A  useful  measure  of  the  decay  of  the  pressure  field  of  the  cascade 
going  away  is  the  ratio  of  the  disturbance  velocities 


(72) 


found  at  distances  x  and  X2  from  the  cascade.  Using  (66)  and  (67), 
this  rado  becomes 


(75)  may  be  simplified  considerably  by  using  a  relation  between  cOqj  and 
Mqq  implied  by  (67).  The  end  result  is 


128 


V.  -  V 

2  m 

V  -  V 

1  « 


■( 


exp  (o)‘2) 
exp  (co^) 


exp  (  to  ’2)  exp  (2^  ) 
exp  (  a>  *x)  exp  (2  co«>7 


(74) 


For 

a)'2« 

0)1 !  «  1, 

a  further  simplification  is  possible: 


v2  ‘  v«  ®  2 


(75) 


v  -  v 

1  « 


CO 


stream;  then,  retaining  the  restriction  (68), 

B  -  l/4  B  (cos  20^  -  i  sin  2^  ) 


d(x  +  iy)  1 _ 

s  “2* 


] 


or 


2*  (x2  -  xx)  /s  = 


(76) 

d  £  and  £ 
from  the  free 

dco  1 

CO  ' 

(77) 

V®  V 

(78) 

[l  -  2  (l/2  exp  o^)2  cos  2^]  ln(o>  '2/co  '  ^ 
1  -  (l/2  exp  co^) 

(1/2  exp  co^)4  ln(  co  '2/co 

s - 


1  -  (l/2  exp  cO 

Inverting  (78)  and  substituting  in  (76), 


I^Jk  =  expifAlV  _ i_- ufe  W _ _ 

V1  y*>  8  l-2(l/2  exp  co^p)2  cos  20^  +  (l/2  exp  o^) 

(79) 

which  is  equivalent  to  equation  (64),  Appendix  II,  for  the  limiting 
case  of  zero  Mach  number. 
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APPENDIX  IV.  ACOUSTIC  PRESSURE  FIELDS  DUE  TO  INTERFERENCE 


BETWEEN  SUCCESSIVE  BLADE  ROWS 


ROTOR  -  STATOR  INTERFERENCE  ANALYSIS 

Several  investigations  (references  3>  15>  18)  Have  studied  the  possi¬ 
bility  of  production  of  propagating  sound  waves  due  to  rotor  -  stator 
interference.  All  three  of  these  analyses  are  quite  similar  and  lead 
to,  among  other  results,  the  concept  of  a  cutoff  frequency  given  by  a 
rather  complicated  expression,  involving  the  radius  ratio  of  an  annular 
configuration  and  one  or  more  Bessel  functions,  usually  of  high  order. 
Some  of  these  analyses  have  been  supported  by  experimental  investiga¬ 
tions  confirming  the  existence  of  a  cutoff  frequency.  The  extent  of 
experimental  confirmation  has  usually  been  rather  slight.  Other 
experiments  often  show  that  rotor  -  stator  combinations  are  noisy  over  a 
wide  speed  range,  with  no  particular  speed  at  which  the  npise  produc¬ 
tion  changes  quantitatively.  Consequently,  a  tentative  conclusion  was 
reached,  that  the  interference  model  may  have  been  unrealistic  or 
that  the  sophistication  of  the  analysis  might  have  hidden  the  essential 
meaning  of  the  phenomena  being  studied. 

A  different  model  is  proposed,  which  will  permit  a  study  of  sound  wave 
propagation  due  to  interferences  between  adjacent  blade  rows,  without 
dependence  on  a  cylindrical  geometry.  Analysis  of  the  model  may  be 
separated  into  two  parts:  (l)  production  and  propagation  of  a  cylin¬ 
drical  wave,  starting  from  a  cylinder  the  minimal  size  of  the  stator 
leading  edge  (this  analysis  is  made  for  a  stator  following  a  rotor, 
but  is  expected  to  apply  also  to  a  rotor  following  an  inlet  guide  vane 
row);  and  (2)  interaction  and  combination  of  cylindrical  wave  fronts, 
with  periodically  spaced  and  time-sequenced  origins,  to  form  plane 
wave  fronts  which  can  propagate  without  attenuation. 

The  wave  front  chosen  here  as  a  fundamental  component  is  cylindrical 
to  represent  the  signal  propagation  from  an  interference  condition  over 
the  entire  blade  span  at  the  same  time.  This  is  inconvenient  in  that 
the  cylindrical  wave  does  not  have  a  simple  acoustic  potential  like 
those  for  plane  and  spherical  waves  (reference  6,  pages  127  “  128, 

163).  However,  the  techniques  of  the  unsteady  flow  method  of  charac¬ 
teristics  (reference  l^,  pages  95-5  et  seq)  may  be  adapted  to  study 
the  propagation  of  a  cylindrical  wave  with  the  initial  conditions:  zero 
particle  velocity  at  zero  time  for  all  radii;  sinusoidally  varying 
particle  velocity  at  the  initial  radius  for  one  period,  followed  by 
zero  particle  velocity  thereafter.  The  FORTRAN  source  list  for  this 
computer  program  follows  the  analysis.  Realistic  numerical  results 
from  this  program  have  not  yet  been  obtained. 

i 

J 
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Figure  61.  Equal  Rotor  and  Stator  Blade  Numbers. 


The  summation  of  cylindrical  wave  fronts  to  form  a  plane  front  requires 
a  determination  of  the  proper  time  relationship  between  reference 
points  on  interference  signals  originating  on  successive  stator  vanes. 
The  simplest  situation  occurs  when  the  number  of  rotor  blades  is  the 
same  as  the  number  of  stator  vanes  (Figure  61)  .  Then  interferences 
take  place  at  the  same  instant  on  all  stator  vanes.  The  symmetry  of 
the  situation  indicates  that  the  individual  signals  will  combine  tc 
form  a  plane  wave  front  parallel  to  the  blade  rows,  which  propagate 
in  the  axial  direction  away  from  the  blade  rows.  The  next  more  complex 
situation  occurs  when  the  number  of  rotor  blades  is  slightly  different 
from  the  number  of  stator  vanes,  and  when  the  time  required  for  a 
signal  to  travel  from  the  leading  edge  of  one  stator  to  the  leading 
edge  of  the  next  stator  is  constant. 


Figure  62.  Varying  Rotor  and  Stator  Blade  Numbers. 
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In  Fignire  62  the  number  of  rotor  blades  is  greater  than  the  number  of 
stator  vanes,  in  the  ratio  5  to  4.  Interferences  occur  later  on  each 
successive  stator  vane  in  the  direction  of  rotation,  by  an  interval 
one-fifth  of  the  time  required  for  one  rotor  blade  to  move  from  one 
stator  vane  to  the  next.  If  this  time  delay  is  less  than  the  travel 
time  of  a  signal  from  one  stator  leading  edge  to  the  next,  the  signal 
continuing  down  a  stator  vane  passage  will  be  a  combination  of  the  two 
signals.  If,  however,  the  first  signal  reaches  the  next  leading  edge 
before  the  signal  starts  there,  there  will  be  no  intersection  between 
the  waves.  Instead  the  signals  will  travel  independently. 

Although  none  of  the  cutoff  frequency  analyses  states  this  explicitly, 
it  seems  likely  that  this  lack  of  wave  intersection  may  be  the  mecha¬ 
nism  behind  the  cutoff  frequency.  Some  qualitative  study  has  suggested 
that  wave  interactions  leading  to  plane  wave  formation  may  occur 
between  signals  originating  farther  apart  in  time.  Quantitative  confir 
mation  of  this  possibility  is  one  of  the  objectives  of  the  computer 
program  that  follows.  Such  confirmation  might  inqjly  a  desirability  of 
revising  the  cutoff  frequency  concept. 


THE  METHOD  OF  CHARACTERISTICS  FOR  CYLINDRICAL  WAVE  PROPAGATION 


Propagation  is  based  on  the  following  equations  (reference  4,  pages  14, 
15,  28,  29)  relating  dependent  variables  (pressure,  p;  density,  p  ; 
sound  velocity,  c;  and  particle  velocity,  u)  in  terms  of  independent 
variables  (time,  t;  and  radius,  r). 

Continuity: 


t-  +  F 

(80) 

Momentum: 

du  ,  du  ,  1  dp  _ 

-rr  +  u  -3—  +  —  -3*-  =  0 
dt  dr  p  dr 

(81) 

Isentropic  flow  of  a  perfect  gas: 

p  =  lip7 

(82) 

Definition  of  sound  velocity: 

2  dp  . ,  7*1 

c  =  =  k7p 

(83) 
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t 


or 


(  7  -  1)  (  y  -  1)  d  lnp 


Using  (83),  (8l)  becomes 

du  ou  2  8 

"5t  +  u  “Sr  +  #! 


dr 


(Inp  )  =  0  . 


(80)  may  also  be  written 


ife  +  -£.  ifcl  +  u  « 0 

dt  r  dr  dr 


or 


dlnp  dlnp  u  du  _ 

■ST— +  u-^rL-  +  F  +  “5?=  0  ' 


(83A) 

(84) 


(85) 


The  potential  velocity  and  density  are  to  be  continuous,  s ingle -vaJued 
functions  of  radius  and  time,  so  that  an  exact  differential  formulation 
is  possible: 


d  lnp 


dr  + 


(86) 


d  u 


du 
=  ~3r 


du 

"3t 


dr  +  -tt  dt 


(87) 


Equations  (84),  (85),  (86),  and  (87)  are  to  be  taken  as  simultaneous 
equations  for  the  partial  derivatives  dlnp  ,  dlnp  ,  du,  du.  A 
typical  derivative,  du,  may  be  expresseS  in  determinant f orm^ a 


as 


“3F 
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-u 

r 


du 

d(lnP  ) 


d  u 

T7 


u 

1 


1 

0 

dt 

0 


1 

0 


u 


c 

u 

0 


0 

1 


0 

1 

0 


principle  of  the  method  of  characteristics  is  that 
uniquely  defined  if  the  denominator  vanishes;  i.e.,  if 


d  u 


(88) 


may  not  be 


dr 

dt  =  uic 


(89) 


For  present  purposes,  attenuation  will  be  confined  to  cases  where  u 
may  be  disregarded  in  comparison  with  c,  so  that 


dr 

dt 


=  +  c 


(90) 


Thus,  lines  along  which 
r  =  r  +  c  t 


(91) 


8  u 

are  lines  across  which  there  may  be  a  discontinuity  in  — yy,  or  any 
other  of  the  partial  derivatives,  or,  put  in  a  different  way,  lines 
along  which  an  acoustic  signal  is  transmitted. 

8  u 

If  u  is  to  be  continuous  along  a  line  such  as  12,  — yy  may  not  be 
infinite,  so  the  numerator  of  (88).  must  also  vanish.  Using  the  plus 
sign  in  (90), 
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0 


-u 

du 

d(lnp  ) 


1  c2  0 

0  u  1 

dt  0  0 

_  dt  , , 

0  c  dt 


=  0 


or  malting  use  of  (83A), 


(92) 


du  + 


7-1 


dc 


uc 

r 


dt  =  0 


(95) 


(93)  defines  the  relation  between  changes  in  c  and  changes  in  u  along 
lines  such  that  radius  increases  with  time.  Alternatively,  if  the 
minus  sign  in  (90)  had  been  used  in  (92),  a  result  would  have  been 

-du  +  ~-j-  dc  +  --  dt  =  0,  (94) 

which  gives  the  relation  between  changes  along  lines  such  that  radius 
decreases  with  increasing  time. 


The  process  of  exploring  the  propagation  of  a  wave  involves  specifica¬ 
tion  of  u  on  all  radii  at  zero  time,  and  over  all  time  at  the  initial 
radius.  Then  (93)  and  (94)  provide  the  necessary  relations  to  construct 
the  field  at  all  times. 
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111*0 
LANE=1 
JLANE*1 

ARG1=6.283185/FlOATF INTIME) 

ARG*ARG1 

JSTOR (l.JFIR ST )*JFIRST+I11 
JST0RI2»JFIRST)=JBEF0R+I11 
JSTOR(3,JFIRST)=0 
JST0RI4.JFIRST )=0 

STCREll,jFIRST)*STOREll.JBEFOR)+OR 
STORE  1 2. JF IRST )=ST0REI2 .JBEFOR )+DR 
ST0REI3.JFIRST )*0. 

STORE (4. JF IRST ) =0. 

J*JF IRST 
Jl*JBEFOR 
J2*0 
J3=0 

JJl* Jl+ JLANE 

5  IF  I JSTOR ( 4*  JJ1  ) ) 9*9*6 

6  J2*J 
J3*J  1 
J*J  +  1 
J1*J1+1 
CALL  CYLWVE 

7  GOTO  5 

9  JSTOR (1«J+1)=J +1+111 
JSTOR (2  «  J  +  l ) =0 
JSTORI3.J+1 J=J1+I 11 
JSTOR (4»J+1)=J+I11 
STORE (1»J+1)=1. 

STORE(2»J+l )=STORE(2  *J) ♦OR 
GOT  0(13,12.12.12) .LANE 

12  STOREI3.J  +  1  )=G. 

GOTO  14 

13  STORE  I  3 .J  +  l )=SINF(ARG) 

14  STORE (4 .J  +  l ) =STORE 1 4  » J ) +GAM1* (  STORE (3»J+1)-ST  ORE  1 3  » J ) ♦ C  STORE  C  3  »  J  + 1 

1  )  /STORE  1 1 . J+l ) +STORE ( 3 . J ) /STORE  1 1 . J ) )*DR/2. ) 

GOTO  (10.16.17,17) .LANE 

10  ARG=ARG+ARG1 
IFIARG-6. 283)11, 15.15 

15  LANE=2 
ARG  *  0. 

11  JBEFOR*JFIRST 
JFIRST=J+2 
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GOTO  4 

16  LANE=3 
JLANE*2 

I F (RATIO! 1 1  — 1 • >39.39.22 
22  LANE=4 
GOTO  11 

17  JSTORI2.J+1)  =  J2  +  1  +111 

I  F(  STORE  (  1. J)-RATIO(KPRINT)  >24.18.18 

39  JJ  =  J+l 
GOTO  40 

18  JJ  *  J 

40  KLIST  =JJ-JF IRST  +  1 

RAT  IO(KPR  I  NT  »*  STORE  ( 1 .  JJ  >  *RA0 
KK=KLIST 

19  JLIST(l»KK)=JSTOR(l»JJ) 

JL  1ST  (2  .tOC>  =  JSTOR(2  .  JJ > 

JLIST ( 3  »KK ) = JSTOR ( 3  » JJ ) 

JLIST ( 4  .  KK ) * JSTOR ( 4 . JJ ) 

ALIST  (  5  .ICICI  =  ST  ORE  ( 1  .JJ) 

ALIST ( 6 .KK ) * STORE ( 2 . JJ  > 

ALIST (7 *KK ) =  AL 1ST ( 6 .KK > -AL I  ST  1 5 .KK > 

ALIST (8 .KK) = STORE ( 3 . JJ ) 

ALIST ( 9.KK) = STORE (4 • JJ) 

AL  I  ST  ( 10  .  KM  =AL  I  ST  (  8  .IOC)  *AL  I  ST  (  5  .  KK  ) 
CTR(KK)=ALIST(7.KK)*RAD/S0VEL 
I F ( KK— 1 >21.21.20 

20  JJJ=JJ 

25  JJJ-JJJ-1 

I F( JSTOR ( l.JJJ)-JSTOR( 3. JJ) >26.27.25 

26  CALL  ERROR 

27  JJ=JJJ 
IOC=ICIC-l 
GOTO  19 

21  WRI TE ( 3. 1CU0 ) SDVEL. GAMMA . RAD . DR. DT.TIME.I (ALIST (K.KK) . 

1  K=l. 10) ,KK=1. KLIST) 

1C  PR  I  NT  =KPRI NT  ♦  1 
DO  41  KK=1. KLIST 

41  BLIST(KK)=ALIST(10.KK) 

CALL  FNWTB«6.1.BLIST,KLIST.IERCOD) 

IF( IERCOD-3 >42.43. 43 
43  WRITEI3.1001 )KPR INT . IERCOD 
SENSE  LIGHT  2 

42  GOTO (11.11.22.23.38) .LANE 
23  IF(KPRINT-NPRINT)24.24.45 

24  IF( JSTORI2.J2) >11.11*8 
8  JJJ=1 

JHOLD= JST  OR ( 1 . 1 ) 

NMOVE=l 

JMOVE=JJJ+NMOVE 

28  DO  29  1=1.4 
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JSTOR ( I  * JJJ ) = JSTOR < I .JMOVE) 

29  STORE ( I  * JJJ  j “STORE ( I .JMOVE ) 

JJJ»JJJ+1 

JMOVErJHOVE+1 

IF ( JMOVE- J) 32 «  32  >33 

32  I F (JSTOR ( 3* JMOVE ) >28*28*30 

30  IF ( JSTOR ( 3. JMOVE l-JHOLD) 28. 31.28 

31  JHOLD=JSTOR ( 1. JMOVE  ) 

JMOVE *JMOVE+l 

GOTO  28 

33  1 1 1* JSTOR ( 1 • JJJ-1 ) - JJJ+1 
JF IRST  =  JJ  J 

34  JJJ=JJJ-1 

IF (JSTOR (4. JJJ I >35.35.34 

35  J8EF0R* JJ J 

DO  36  I = JFIRST .J 
DO  36  11=1.4 
STORE  1  1 1 • I ) =0. 

36  JSTOR (1 1. 11=0 
GOTO  4 

37  LANE=5 
GOTO  18 

38  CALL  EXIT 

45  CALL  FNWTBI6.2.CTR.KLIST. RATIO. NPRINT.IERCOD) 

IF( IE RC 00-3)46. 47. 47 
47  WRITE ( 3  *1001 JKPR I NT 
GOTO  1 

46  IFtLTAPE 1 1.80.80 

80  REWIND  6 
KERR0R=2 

DO  79  1=1.7320 
79  STORE ( I ) =0. 

DO  81  KPRINT=i.NPRINT 

CALL  FNRTB( 6 .1 .U( l.KPRINT ) .KL I  ST . IERCOD ) 

IF( IERCOD -6>81.82. 82 

82  WRITE(3. 1002 IKPRINT . IERCOD 
GOTO  1 

81  CONTINUE 

CALL  FNRTB(6.2.CTR.<L I  ST. RAT  10. NPRINT.IERCOD) 

IF( IERCOD-6 184.83  .83 

83  WRITE(3»10U2 1XPRINT . IERCOD 
GOTO  1 

*  SPACE  IS  THE  DISTANCE  BETWEEN  TWO  STATOR  VANES. 

84  BLADR=NBL ADR 
BLADS=NBLADS 

NDIST=LCOUNT (1 .100.DIST.1 ) 

*  TMEVNT  REFERS  TO  THE  TIME  BETWEEN  INTERFERENCE  EVENTS  ON  ADJACENT 

*  STATOR  VANES. 

*  CYCLE  REFERS  TO  THE  TIME  BETWEEN  PASSAGE  OF  SUCCESSIVE  ROTOR 

*  BLADES  PAST  A  POINT. 
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*  TIME  REFERS  TO  PASSAGE  OF  A  BLADE  WAKE  PAST  A  POINT. 

TMEVNT« ( BLAr R/BLADS-1 . ) *SPACE/BLDSPD 
CYCLE«SPACh*BLADS/BLADR/BLDSPD 
TIME*WAKE/BLDSPD 
DO  147  KriST«l»NDlST 
NAXIS*DIST(KDIST  >/SPACE*20. 

NAX I S*XMAXOF (NAX I S  »  30 ) 

*  CYCTIM  GIVES  THE  TIME  AFTER  THE  START  OF  SOME  SINE  WAVE  SIGNAL 

*  AT  THE  SOURCE  WHEN  A  SIGNAL  REACHES  THE  POINT  BEING  STUDIED. 

*  CYCTMM  GIVES  THE  TIME  WHEN  THE  END  OF  THE  SIGNAL  REACHES  THE  POINT. 

C YCT IM=M0DF <  DI ST  <  KD I  ST ) /SDVEL .CYCLE  ) 

LANE  1*1 

CYCTMM*CYCT IM+T IME 
'F(CYCTMM-CYCLE 1 101. 101.100 

100  CYCTMM*CYCTMM-CYCLE 
LANE1*2 

101  C YCTX*0, 

DO  108  JCYC=1 » 100 
VY ( JCYC ) =0. 

GOTO! 105.102) . LANEl 

102  I F ( CYCTX-CYCTMM) 103.103.104 

103  CALL  BI LAG( CYCTX-CYCT IM+CYCLE  »D 1ST ( KOI  ST )*SPACE .UCOMP ) 

VX( JCYC)=UCOMP*RAD/DIST(KOIST) 

GOTO  108 

104  LANEl *1 

106  VX( JCYC )*0. 

GOTO  108 

105  I F (CYCTX-CYCTI M) 106. 106 .107 

107  LANE1*2 
GOTO  103 

108  CYCTX*CYCTX+CYCLE/100. 

150  DUM  *  DUM 

DO  133  KAXI S* 1 .NAX I S 
AXIS*KAX I S 

D I STT  *SORTF I DI ST  I KD I  ST ) **2+AXI S**2 I *SPACE 
CYCT IM=MODF ( DI STT/SDvEL .CYCLE ) 

SINANG*SPACE*AXIS/DlSTT 
COSANG=DIST (KDISTI/DISTT  *  SPACE 
T IMPLS*CYCT IM+TMEVNT*AXIS 
T IMMNS*CYCT IM-TMEVNT  *AX IS 
148  I F ( T I MPLS  1109.110,110 

109  T IMPLS*T IMPLS+CYCLE 
GO  TO  148 

110  IF(TIMPLS-CYCLE)112.111.111 

111  T IMPLS*T IMPLS-CYCLE 
GOTO  HO 

112  IF ( T IMMNS  >113.114, 114 

113  T IMMNS*T I MMNS+CYCLE 
GO  TO  112 

114  I F ( T IMMNS-CYCLE 1116. 115. 115 


141 


115  T IMMNS=T I MMNS-C YCLE 
GOTO  114 

116  CYCTX*0. 

T 1MPLX=T 1MPLS+T 1ME 
LANE l*i 

IF (T1HPLX-CYCLE 1118.118.11 7 

117  LANE1*2 

T IMPLX*T 1 MPLX-CYCLE 

118  DO  124  JCYC-1.100 
GO TOI122.119) .LAN El 

119  IF(CYCTX-T IMPLX  >120.120.121 

120  CALL  BILAG(CYCTX-TIMPLS«CYCLE»DISTT.UC0MP) 

VX  ( JC  YC  )  *  VX  (  JC  YC  )  4-UC0MP*RAD/D  I  ST  <  KD I  ST )  *COS ANG 
VY ( JCYC ) * VY ( JC  YC  ) -UCOMP*RAD/DI ST (HOIST) *S:NANG 
GOTO  124 

121  LANE1*1 
GOTO  124 

122  I F( CYCTX-T IMPLS ) 124.124.123 

123  LANE1-2 
GOTO  120 

124  CYCTX*CYC7X+C YCLE/ 100. 

CYCTX«0. 

TIMHNX*TIHHNS+TIME 

LANE1=1 

IF(TIMMNX-CYCLE) 126.126.125 

125  LANE  1*2 

T IMMNX=T IMMNX-CYCLE 

126  DO  132  JC.C=1.100 
GOTO (130.127) .LANE1 

127  IF(CYCTX-T IMMNX (128.128.129 

128  CALL  BILAGICYCTX-TIMHNS+CYCLE.DISTT .UCOMP) 

VX  (  JCYC  )  =  VX  (  JC  YC  )  +UCOMP*RAD/DI  ST  ( ICD  I  ST  )  *COSANG 
V  ( JCYC  »  =VY  (JCYC)  +(JCOMP*RAD/DI  ST  (  KD  I  ST  )  *SI  NANG 
GOTO  132 

129  LANE1=1 
GOTO  132 

130  I F (CYCTX-T I MMNS ) 132.132.131 

131  LANE  1=2 
GOTO  128 

132  CYCTX=CYCTX+CYCLE/100. 

133  CONTINUE 
VRMAX*0. 

DO  140  JC YC=1. 100 

VR( JCYC ) =SQRTF ( VX ( JCYC  » **2+VY ( JCYC ) **2 ) 
IF(VXIJCYC) >134,135.136 

134  VANGI JCYC 1=3.14159 
GOTO  137 

135  VANGI JCYC)* SIGNF(VY( JCYC) .1.5708) 

GOTO  138 

136  VANGI JCYC )=0. 
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137  VANG(JCYC)»VANG( JCYC)+ATANF(VY( JCYC  )/VX(JCYC) ) 

138  VRMAX*MAX1F(VRMAX,ABSF (VRIJCYC ) )) 

IF ( VRMAX-ABSF ( VR ( JCYC ) ) 1 139. 139 , 140 

139  LVRMAX* JCYC 

140  CONTINUE 

DO  141  JCYC=1» 100 
X( JCYC+1 I*.01*FLOATF( JCYC) 

VNORM ( JC  YC  )  *  VR  ( JC  YC  )  »COSF  ( V  ANG  (  JC  YC  ) -VAN6 ( LVRMAX ) J 

141  VTANGI JC YC I *VR (JCYC ) »S1NF ( VANG ( JC YC ) -VANG ( LVRMAX ) ) 

XII)  *  0. 

VX(101)«VX 
VY ( 101 ) *V Y 
VR(10l)«VR 
V  ANG  ( 1 0 1 )  * V  ANG 
VNORM ( 101 )= VNORM 

142  VTAN6I101 )=VTANG 
CALL  ORG(0.,5.  ) 

CALL  SCALEMI 1. 6.0 . .+VRMAX.-VRMAX. 16 . »10. I 
CALL  XAXISAIO. *0. .BCD! 1 ) .-4. 16. .0. ) 

CALL  YAX1SAI0. .-5 . .BCD! 2 ) .57 , 10 . *90. ) 

CALL  SYMBLAI .75,10. ..21, TITLE »0.,46 I 
CALL  SYMBLAI 10. ,9. 5 . . 14.LABEL1 . 0. .42 ) 

143  XPLOT*BLDSPD/SDVEL 

CALL  NUMBRA ( 15.2  *9.5.. 14 , XPLOT • 0. .3  ) 

CALL  SYMBLAI 10. . 9.2 • . 14.LABEL2 . 0. .42 ) 

144  XPLOT-WAXE/SPACE 

CALL  NUMBRA 1 15.2 .9.2 . . 14. XPLOT .0. .3 ) 

CALL  SYMBLA (10. .8. 9*. 14.LABEL3 . 0. .42 ) 

CALL  NUMBRA! 15.2 . 8. 9 . . 14 .BLADR . 0. .- 1 ) 

CALL  SYMBLA (10. .8. 6. . 14 .LABEL4.0. .42 ) 

CALL  NUMBRA! 15.2 .8.6..14.BLA0S.0* .-1) 

CALL  SYMBLAI 10. .8.3. . 14.LABEL5 .0. .42 ) 

145  XPLOT«DIST(KDIST J/SPACE 

CALL  NUMBRA ( 15 .2 . 8. 3 . . 14 .XPLOT , 0. .2 ) 

CALL  SYMBLA (10. .8... 14.LABEL6.0. .42 ) 

146  XPLOT *VANG( LVRMAX ) *57.296 

CALL  NUMBRA! 15.2 .8.0 .. 14, XPLOT ,0. .2  I 
CALL  LINEIX. VNORM. 101,1) 

CALL  LINEIX, VTANG, 101. 1) 

CALL  SYMBLAI X( LVRMAX) -.25. VNORM (LVRMAX )+. 05*. 14.LABEL7I1 ) , 

1  0. ,5 ) 

CALL  SYMBLAI X( LVRMAX )-. 25 .VTANG (LVRMAX )+.05 , . 14. LABEL 7 ( 2 ) , 
1  C.,5) 

GOTO  (147,147,38,38) .KERROR 

147  CONTINUE 
GOTO  1 

1000  FORMAT 
RESTORE 

CYLINDRICAL  ACOUSTIC  WAVF.  PROPAGATION  ANALYSIS 
SPACE 
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SOUND  VELOCITY 

GAMMA 

INITIAL  RADIUS 

-0PF1 

-F2 

-1PG3 

FT 

RADIUS  INCREMENT 

-1PG3 

(DELTA  Rl/R 

TIME  INCREMENT 

-G3 

SEC. 

TOTAL  WAVE  OURATION 

-G3 

SEC 

SPACE 

PN  PNRC  PNDG  PNLC 

R/RO 

TIME  (CT-R) 

U  (C-CO)  RO/R 

X»U 

SPACE 

-X  -I  -I  -I  -I 

-0PF4 

1 

■n 

* 

I 

■n 

* 

1 

“n 

* 

I 

"n 

** 

XF4 

REPEAT  1 
END  OF  FORMAT 

1001  FORMAT 

TAPE  RECORD  NO  -I  WAS  NOT  WRITTEN  SUCCESSFULLY »  CODE  -I 
END  OF  FORMAT 

1002  FORMAT 

TAPE  RECORD  NO  -I  WAS  NOT  READ  SUCCESSFULLY »  CODE  -I 

END  OF  FORMAT 

END 
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C  erLWVE  CHARACTERISTIC  lNTERSEC-'  .‘ON  SUBROUTINE 
SUBROUTINE  CYLWVE 

DIMENSION  U ( 60  *1 00 )  »CTR 160 1 >01 ST  1 100) .STORE  14.1830) .JSTOR! 4. 1830 ) 

1  .RATIO! 100) 

COMMON  U. JS TOR. CTR. DIST. BLOSPD. SPACE. NBLADR .NttLADS. WAVE. RAD. KPRI NT 

1  .NPRINT .NDIST .KDIST .RATIO. J.Jl .J2.J3.Ill »0K .bAMMA .oAMl.SOVcL 

2  .GAM11 
EQUIVALENCE  (U.STORE) 

•  J1  REFERS  TO  L  CHARACTERISTIC 

•  J2  REFERS  TO  R  CHARACTERISTIC 

•  J3  REFERS  TO  DIAGONAL  POINT 

•  J4  REFERS  TO  NEXT  L  CHARACTERISTIC  POINT 

•  J5  REFERS  TO  NEXT  R  CHARACTERISTIC  POINT 

JSTOR! 1 . J )*J+I 11 
JSTOR  (2  .  J  ) * Jl* I 1 1 
NCOUNT  « 1 

JSTOR ( 4. J )*J2+Ill 
IF( J3I1.1 .3 
1  JSTORI3.J )*0 

STORE ( 3 . J )  *  I  STORE  I  3 » J1 ) ♦STORE ( 3 • J2 ) ) /2 • 

STORE (4. J )  >  ( STORE !4.J1)+ST0RE(4.J2) ) /2 • 

GOTO  5 

3  JSTOR (3 .J ) * J3+I 11 

STORE  ( 3 . J )  >  STORE ( 3. J1 l+STORE ( 3 .J2 1 -STORE ( 3  * J3 ) 

STORE ( 4* J )*STORE ( 4. J1 ) *STORE (4. J2 )-ST0RE(4, J3  ) 

6  STORE (l.J )*STORE( 1.J3) 

STORE  1 2 .J )* STORE ( 2  *  J  3 ) +2  »*DR 
IF(JST0R(2.J1) 15.5.11 

5  C1*GAM1 1#DR* (STORE ( 3. J) /STORE! 1 . J ) +STORE ( 3 . J1 ) /STORE! 1 .J1 )) /2. 

1  ♦GAM1*ST0RE( 3 . J1 l+STORE ( 4* J1 ) 

GO  TO  13 

11  J4* JSTOR ( 2. J1 ) -I 1 1 
IF(J4)5»5.12 

12  C1*GAM11*DR» I «416667»STORE( 3 » J I /STORE ( 1 »J)+ .66666 7*STORE ( 3 . J 1 ) 

1  /STORE! l.Jll-.083333»STORE( 3. J4) /STORE (1.J4) ) +GAM1 ‘STORE ( 3 . Ji ) 

2  ♦ST0REI4.J1 ) 

13  IFIJSTORI4.J2) )14. 14*16 

14  C2*GAM11*DR» (STORE! 3 • J I /STORE ( 1 • J )♦ STORE (3  »  J2 ) /STORE ! 1 . J2 ) ) /2. 

1  ♦GAM1 l»STOREI 3 • J2 l+STORE (4 . j2 ) 

GO  TO  15 

16  J5«JST0R(4.J2)-I11 
IF  ( J5H4.14.17 

17  C2*GAM11*DR» ( «416667*STORE « 3 .J ) /STORE! 1 * J)+ .66666 7*STQRE ! 3 . J2  ) 

1  /STORE! 1 »J2)-.083333*STORE! 3. J5  J/STORE!  1 . J5) J+GAM1 1 

2  »STORE(3»J2)+STORE(4.J2) 

15  CALL  DETER(C1.C2»GAM1»GAM11»1..1..CR*UR) 

7  IF(ABSF» CR-STORE! 4. J ) H-ABSF ( UR-STORE!  3  •  J ) l-l.E-5  19.9.0 

8  ST0REI4.J )«CR 
STOREI3.J )«UR 
NCOUNT «NC0UNT+1 
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IFINCOUNT-15 )5  »5» 10 
10  CALL  ERRORA 
9  RETURN 
ENO 

CB1LAG  BIVARIATE  LINEAR  INTERPOLATION.  CYLINDRICAL  WAVE  INTERFERENCE 
SUBROUTINE  BILAGIX.Y.Z) 

DIMENSION  UI60.100J .CTRI60) »DIST( 100) .STORE (4*1830) .JST0RI4* 1830 1 
1  .RAT IO( 100 ) 

COMMON  U • JSTOR «CTR. D 1ST. BLDSPD. SPACE. NBLADR.NBLADS.WAVE. RAD. ICPRINT 
1  .NPRINT.NDIST.ICDIST  .RAT  10.  J.  J1  . J2.J3.  II 1  »DR. GAMMA .GAMl.SDVtL 

EQUIVALENCE  (U. STORE) 

1*1 

J*1 

1  1 F ( X-CTR ( 1+1 ) )4«4.2 

2  1*1+1 

1 F ( 1-59 11.3.3 

3  CALL  ERROR 

4  1F(Y-RATI0< J+1II7.7.5 

5  J*J+1 

IF ( J-NPR1NT 14.6,6 

6  Z  =  UU.JI  +  IX-CTRm)/ICTR«l+ll-CTR«m*(U(l  +  l,J)-UII.J) 

GOTO  8 

7  DEN*  (CTR(  1  +  H-CTRI  1  )  *  *  C  D1ST  C  J+l  )-DIST(  J)| 

CLB*V-01ST( J| 

CUB«01ST(  J+D-Y 
CLH*X-CTR(H 
CRH*CTR ( 1 +1 l-X 

Z*(CUB»CRH#U < 1 »J) +CUB»CLH*UI 1  +  1 . J  J +CLB*CRH*U( I  * J+l ) 

1  +CLB*CLH*U ( 1 +1 » J+l ) l/DEN 

8  RETURN 
£N0 
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ANALYTICAL  INVESTIGATION  OF 


RADIATION  CHARACTHUSTICS 


FROM  A  FAN  FACE 


The  approach  used  is  that  of  a  piston  with  a  varying  particle  velocity. 
The  particle  velocity  was  assumed  to  be  a  function  of  the  radial  posi¬ 
tion  and  the  difference  in  the  number  of  rotor  and  stator  blades.  If 
the  number  of  rotor  and  stator  blades  is  the  same,  the  source  is  essen¬ 
tially  a  piston  radiator,  as  all  the  rotor  -  stator  excitations  are  in 
phase.  When  the  number  of  rotor  blades  exceeds  or  is  less  than  the 
number  of  stator  blades,  the  pressure  fluctuations  are  out  of  phase 
and  the  directivity  pattern  varies  with  an  increase  in  the  difference 
of  rotor  and  stator  blades,  the  wave  length,  frequency,  and  rotor  dimen¬ 
sions.  The  basic  equation  for  the  derivation  was  the  acoustic  pressure 
at  a  point  resulting  from  an  excitation  in  a  plane  a  distance  r '  away 
(reference  6). 

Figure  63  is  a  schematic  of  the  coordinate  system  used  for  the  deriva¬ 
tion. 

dp  -  u3  eX-Ve)*  (95) 

where 

ds  =  element  of  rotor  surface  area, 
dp  3  acoustic  pressure  at  point  P. 

P  =  density  of  medium. 

=  frequency  of  excitation  (  “>  =  x  RPS). 

uj  -  particle  velocity, 

t  =  time. 

c  =  speed  of  sound  in  medium. 

N  3  number  of  rotor  blades, 

r 

BPS  3  revolutions  per  second  of  rotor. 

The  particle  velocity  may  be  assumed  to  have  the  form 

u5  =  UQ  (r)  e  Jn[P  (96) 

where 

n  3  difference  in  the  number  of  rotor  and  stator  blades. 
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Substituting  (96)  into  (95)# 


dp 


Jl>“Uo(r)eJ  n\J/+a>  (t 
2  x  r ' 


(97) 


with 


ds  -  rdrdlj/  . 

From  Figure  63,  the  following  may  be  determined: 


r'2  =  (r  coslj/)2  +  ( R  sin  6  -  r  sln\j/ )2  +  (r  cos  6  )2 
=  r2  +  R2  -  2rR  sin  0  sin \j/ 

2 

r'  =  R(1  +  -  2  |  sin  a  sin 


«nd  If  (|)2<SC  | 


r' 

=  R(1 

-2|sln 

0  sin  if/)1'2 

or 

r' 

R(1 

-  ^  sin  6 

sin  \J/) 

r' 

(R  - 

-  r  sin  0  sin  \j/ ) . 

Substituting 

equation  (98)  into  (97)# 

dp^Jpw  UQ 
2  x  R 

eJ 

(t  -  R/c)t 

J  (n\J/  +  r 

c 

(98) 


(99) 


Note  that  the  last  term  of  equation  (98)  was  omitted  in  the  denominator 
of  equation  (99)#  but  retained  in  the  exponential  term.  This  was  done 
by  assuming  that  the  term 

Jpo>  U0 

2  jt  r ' 

oscillates  about  the  value  1,  which  enables  the  term  r  sin  6  sin  \J/ to 
be  assumed  negligible,  while  the  exponential  terms  oscillate  about  a 
value  of  zero,  which  does  not  allow  the  term  r  sin  0  sin  \l/to  be 
assumed  negligible.  ' 
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From  '^Tables  of  Functions"  by  Jahnke  and  Emde,  equation  (99)  becomes 

P  _  e  j  u)  (t  -  R/c)  J  Uq  Jn  (_cd£  sin  Q  ^rdr  (10o) 

rv  C 

from  which 

u  =  cos  o>  (t  -  R/c)  (  U  (r)  J  (-^-£  sin  0  )  rdr 

Rc  J  o  n  c  (101) 

I  =  ^ —  |fu  (r)  J  (oj  r  sin  0  )  rdr  2.  (102) 

4R2c  [J°  n  T  J 

Equations  100,  101,  and  102  were  programmed  to  provide  an  output  of 
sound  pressure  level  relative  to  an  assumed  reference  pressure  level  as 
shown  in  the  following  FORTRAN  source  list. 
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FORTRAN  SOURCE  LIST  FOR  RADIATION  FROM  A  BOTCH  ANNULUS 


*FANAma§Y8R  imtl IS»^W?i«i5K88ECmYi »  .SOPARI  19)  *TritTA(2Q) 

1  .SDPRESI 19 ! 

TABLE  8115(0-377777777777) 

IF  I  SENSE  LIGHT  1)  6.8 

8  RHO  -  .076475 
SDVEL  »  1)16.4 
DO  2  I  •  1.38 

2  RADI  I )  »  BITS 

READ  DIP  RPM .BLADES. RAD »U0» SDVEL. RHO .BEATS 
IE  (  SENSE  LIGHT  1)  8.1 
1  ALBDA  «  SDVEL*9.54929/BLADES/RPM 
NRAD  *  LC0UNTd.19.RAD. 2) 

NR ADI  *  ( NRAD+1 ) /2 

SDREF  *  RHO*SDVEL/4.*RADI2,NRADl)**2/32. 174/685. 399 
REF  «  RAD (2 .NR ADI ) 

PREF  «  RHO*SDVEL/32.l74*REF/2.08848E-3 

DO  3  I  *  l.NRAD 

RADII, I)  *  RADII .1 l/ALBDA 

3  RADI 2 , 1 )  *  RAD (2  *  I ) /REF 
THETA! 1 )  «  0. 

SDPAR(l)  »  0. 

9  NBEAT  -  BEATS 
NBLAD  *  BLADES 
DO  7  I  «  1,19 

THETA!  Id  )  *  THETA!  I-1J+5. 

THET  *  SINFITHETAI I )*. 0174532) 

DO  4  J  a  l.NRAD 

CALL  BESSL 1 1 .RAD  1 1  *  J ) *THLT .NBEAT .NBEAT +5 .BESSEL  1 1 . J ) ) 

IFISENSE  LIGHT  1  )  6.4 

4  CONTINUE 

TRY1  a  RADI  1*1 )*RAD( 2 • 1 ) *BESSEL (1.11 
SDPRESII)  *  0. 

DO  5  J  *  3, N RAD, 2 

TRY2  -  RADI  1 • J-l ) *RADI 2  »J-1 ) *BESSEL 1 1.  J-l ) 

TRY3  a  RADII. J ) *RAD I  2  .  J ) *BESSEL 1 1 . J  ) 

SDPRESI I ) “SDPRES I  I l+IRADI l.Jl-RAOl 1 .J-2 ) )•( .166667*1 TRY1+TRY3)* 

1  .666667*TRY2 ) 

5  TRY1  a  TRY3 
SDPAR(I)  -SDPRES I  I ) **2 
GOTO  7 

6  WRITE  I  3. 1000) THETA! I ) .J. RADI  1 .J ). THET .NBLAD 

7  CONTINUE 

WRITE! 3 .1001 )NBEAT.NBLAD.PREF, RPM .SDREF, ALBDA. RHO. REF, SDVEL. 

1  (RADI l.J),RAD(2*J) .THETA! J) .SDPARIJ) .SDPRESI J) »J«1, 19) 

GO  TO  8 
1000  FORMAT 
SPACE 

BESSEL  FUNCTION  SUBROUTINE  ENCOUNTERED  AN  ERROR  AT  THETA  «  -F3.  J«  - 
X  RADI  1 . J) »  -E3.  SINTHETA  «  -F3.  NBLADES  «  -I 
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END  OF  FORMAT 
1001  FORMAT 
RESTORE 

SOUMO  RADIATION  FROM  A  ROTOR  BLADE  ANNULUS 
SPACE 

NUMBER  OF  BEATS  ■  -I 

NUMBER  OF  BLAOES  >  -1 

XESSURE  «  -1PG3  DYNES/CM2 

RPM  ■  -OPFO 

XTENSITY  «  -1PG3  WATTS/CM2 

REFERENCE  LENGTH  *  -0PF3  FT  (2*PI»MAVE  LENGTH i 

X  «  -FS  LBM/FT3 

REFERENCE  VELOCITY  ■  -F3  FT /SEC 


XTY  ■ 

SPACE 

-FI  FT/SEC 

RADIUS 

PEAK  VELOCITY 

ANGLE  FROM  AXIS 

xntensitv 

SOUNO  PRESSURE 

REFERENCE  VELOCITY 

(DEGREES) 

XR/LBOA ISO 

RAT I0*( R/L8DA) 

-OPFJ 

-F3 

-FO 

X-1PG3 

-63 

REPEAT 

1 

ENO  OF 
END 

FORMAT 

REFERENCE  PR 
REFERENCE  IN 
AIR  DENSITY 
SOUND  VtLOCI 

SOUND  I 
RATIO* I 
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n  n  n 


ceCSSLe2Ss£L  FUNCTION  JIX)  ANO  I « X ) 

BESSEL  SUBROUTINE 
CHANOLER  MARCH  -  TIS  64TIP5 

STEGUN  ANO  ABRAMQWITZ.  MTAC  60  P.252  -  OCT  196? 
SUBROUTINE  8ESSLI ITYPE.X.NMIN.NMAX.BESJI ) 

C  ITYPE  «  1  FOR  UIX) 

C  ITYPE  •  2  FOR  I(X) 

C  STATEMENTS  WITH  100  SERIES  EFN  LUMBERS  00  NOT  APPLY  TO  I(X) 
3  DIMENSION  BESJK10) 

10  XX  '  X 
MAXN-NMAX 

M1NN-NMIN 
NOEL «MAXN -MINN 
IF (NOEL  >11*14*14 

11  NTMP«MAXN 
MAXN=MINN 
MINN«NTMP 
NOEL “-NOEL 

C  NOEL  IS  NUMBER  OF  ORDERS  OESIRED 

14  NDEL»NDEL*1 
IFlXXI43.4l.43 

41  00  42  1-1 .NOEL 

42  BESJUll-O. 

GO  TO  (44.40). ITYPE 

44  IF (NMIN >46*45 .40 

45  BESJ1 1 1 )*1. 

GO  TO  40 

46  I F (NMAX ) 40.47.46 

47  BESJ1 (NOEL)bI. 

GO  TO  40 

48  N=1-NMI N 
BESJl (N)»l» 

GO  TO  40 

43  NOEL 1 -NOEL 
G0T0I15. 116). ITYPE 

15  Cl  •  -1. 

C2  •  2. 

GOTO  17 

116  Cl  »  1. 

C2  •  1. 

17  K  «  ( XMAXOF ( MAXN.XF IXF ( 1«5*XX+. 51 ) ) /2 ) »2+12 
FJP3  •  0. 

FJP2  •  l.E-11 

C  SUMFJ  ACCUMULATES  TENTATIVE  BESSEL  FUNCTION  VALUES 
C  TO  GET  THE  NORMALIZING  FACTOR 
SUMFJ  «  0. 

TWOX  •  2. /XX 
P  »  K*1 

18  FJP1  ■  TW0X*(P+1« )*F JP2*C1*F JP3 
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ABFCN1  «  ABSF(FjPl) 
i F(ABFCNl  -  1.E36 >20*19*19 

19  SENSE  LIGHT  1 
GOTO  40 

20  FJPO  «  TW0X*P*FjPl*Cl*FjP2 
ABFCNO-  A8SF (FJPO ) 

IF(ABFCN0-1.E36)21.19.19 

21  IF  (ABFCN0+ABFCN1-1.E-30 122*23*23 

22  FJP0«0. 

FJP1  «  0. 

23  GOTO (25 *24) *  I TYPE 

24  SUMF J-SUMF J+FJP1 

25  SUMF J*SUMFJ+C2#F JPO 
IF (NOEL) 32 *32*51 

51  I F ( K-MAXN ) 50  *29*  32 

50  IF  (K+l-MAXN) 32  *21*32 

C  THE  FUNDAMENTAL  RECURRENCE  FORMULA 
C  FOR  J.  JP(X)  «  2»(P*1>/X*JP1(X>-JP2(X) 

C  FOR  I*  IP(X)  ■  2»<P*1>/X*IP1(X>-HP2(X) 

27  BESJI (NDEL)-FjPl 
NDEL-NDLL-1 
MAXN«MAXN-1 
IF (X-NMIN 132*29*29 
29  BESJI (NOEL) -FJPO 
NDEL-NDEL-1 
MAXN-MAXN-1 
GOTO  32 

C  THE  SUM  IN  THE  NORMALIZATION  CONSTANT  FOR  J  CONTAINS 
C  ONLY  EVEN  X  IN  JX 

32  IF(X  >34*34*33 

33  P-P-2. 

X-X-2 
FJP3-FJP1 
FJP2-FJP0 
GOTO  18 

34  SUMF J-SUMFJ-FJPO 
GOTO) 36* 135 ) » I  TYPE 

135  SUMF J-SUMFJ/EXPF ( XX ) 

36  DO  37  I-l *NDEL1 

37  BESJ I ( I ) -BESJ I ( I ) /SUMF J 
40  RETURN 

END 
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APPENDIX  VI. 


ABSORPTIVE  SUPFRESSOR  DESIGN  CALCULATIONS 


DESIGN  PROCEDURE  FOR  ABSORPTIVE  LOUVERS 


A.  Determine  percent  open  area,  number  of  louvers,  axial  length  of 
louvers,  and  the  frequency  of  the  noise  source. 

B.  Calculate  ly  (louver  spacing). 

C.  Calculate  X  (wave  length). 


D.  Consult  reference  9>  and  read  the  db  attenuation  per  length  equal 
to  ly,  Aly,  for  the  various  values  of  ly/ ^  for  the  design  percent 
open  area. 

E.  Multiply  Alv  by  — —  to  determine  predicted  noise  reduction,  Mb. 

y 

F.  Read  F  (flow  resistance  parameter)  corresponding  to  design  percent 

open  area.  R, 


Calculate  R^  and,  with  reference  9/  select  proper  absorptive 
material. 


SAMPLE  CALCULATION 


80 
10 

1  ft. 

1070  cps  and  harmonics 

B.  1  =  2.04"  (calculated) 

y 


A.  Assume:  percent  open  area  = 

No.  of  louvers  =■ 

Axial  length  of  louvers  = 
Frequency  of  Source; 
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c. 


f  (cps)  X(ft) 

1070  1.05 

2140  .527 

3210  . 352 

4200  .264 

5350  .210 

D-  VX  “y 

.162  .3 

.322  1.2 

.483  1.9 

.644  2.6 

.81  3 

E.  A1  A(A  db) 

.3  1.76 

1.2  7.05 

1.9  11 

2.6  15.3 

3  17.6 


F.  F  =  .75 


t  =  l/4"  =  .00635  meters  pC  =  406  Mks  Rayles  for  air 

R^  =  4.8  x  10^  Mks  Rayles/m 


From  reference  9, 


select  proper  material  for  this  R 


r 
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